FORMAL DESINGULARIZATION OF SURFACES 
- THE JUNG METHOD REVISITED - 
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' Abstract. In this paper we propose the concept of formal desingularizations 

Ph ' as a substitute for the resolution of algebraic varieties. Though a usual res- 

olution of algebraic varieties provides more information on the structure of 
singularities there is evidence that the weaker concept is enough for many 
1/^ ' computational purposes. We give a detailed study of the Jung method and 

show how it facilitates an efficient computation of formal desingularizations 
for projective surfaces over a field of characteristic zero, not necessarily al- 
' . ' ■ gebraically closed. The paper includes a generalization of Duval's Theorem 

L , _ , J ' on rational Puiseux parametrizations to the multivariate case and a detailed 

' description of a system for multivariate algebraic power series computations. 
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1. Introduction 

Smooth varieties are (in general) well-understood. By contrast (or simply be- 
cause of that) the objects of interest are often singular varieties. From the theoreti- 
cal point of view, a remedy for this situation is the celebrated Theorem of Hironaka 
[20] (or [19] for a more modern treatment) on the resolution of singularities which 
is ubiquitous in algebraic geometry: If A is a variety over a field of characteristic 



Date: February 2, 2008. 

2000 Mathematics Subject Classification. Primary: 14Q10, 14J17, 13J05. 

Key words and phrases, resolution of singularities, algebraic power series, quasi-ordinary 
polynomials. 

The author was supported by the FWF (Austrian Science Fund) in the frame of the research 
project SFB F1303 (part of the Special Research Program "Numerical and Symbolic Scientific 
Computing"). 

1 



2 



TOBIAS BECK 



zero, then there always exists a smooth variety Y and a proper birational mor- 
phism TT : Y X. So for proving theorems and defining birational invariants, 
one can often argue on Y rather than on X and finally transfer the result back 
to the singular variety. This theorem has been made constructive by Villamayor 
[37], Bierstone-Milman [5] and others. There arc also two implementations of the 
dcsingularization algorithm in Singular [16], one by Bodnar and Schicho [6] and 
another one by Friibis-Kriiger [14]. In principal this makes many theoretical results 
algorithmic, but any algorithm relying on a dcsingularization suffers from the high 
computational complexity of the dcsingularization process. There are also special- 
ized constructive methods for the surface case (for a list see [9] and [24, Clip. 2]), 
in particular, the Method of Jung which originates in [22] and has been further 
developed in [38] and [21]. But - to the best of our knowledge - there are no 
implementations available. 

From the computational point of view, it is not always necessary to describe a 
dcsingularization completely. In the case of algebraic curves over a field E, series 
expansions have proven to be an important algorithmic tool. Here the preimage of 
the singular locus w.r.t. a dcsingularization is a finite set of points. The idea is to 
describe the dcsingularization by power series expansions that determine "formal 
neighborhoods" of these points. If E has characteristic zero, Puiseux expansions 
can be used, and the Newton-Puiseux algorithm is implemented in many computer 
algebra systems including Magma [7], Maple and Singular. The latter system also 
has an implementation of Hamburger-Noether expansions [8] that provides a similar 
tool for positive characteristic. Applications are for example the computation of an 
integral basis of the hmction field [36] or Riemann-Roch spaces of divisors [17]. 

The purpose of this paper is to provide a similar tool for hypersurfaces of Pj|, 
where E is a field of characteristic zero. We emphasize algorithmic aspects and 
proceed as follows: In Section 2 wc define formal dcsingularizations for schemes 
of arbitrary dimension. They can be interpreted as sufficiently large sets of local 
paramctrizations by formal maps. Formal dcsingularization offer a lot of flexibil- 
ity because during computation one can always to switch to formally isomorphic 
schemes. Then in Section 3 we show how to compute them for surfaces using the 
method of Jung that depends crucially on the Theorem of Jung-Abhyankar. We de- 
fine and use rational Puiseux paramctrizations whose existence and computability 
we assume for that moment. We also give a description of the algorithm in mathe- 
matical pseudo-code. Finally Section 4 shows in detail how to implement a system 
that represents and computes with multivariate algebraic power series in Magma. 
(Locally smooth systems, as proposed in [1], were insufficient from the complexity 
point of view.) Folklore knew that the concept of rational Puiseux paramctriza- 
tions introduced by Duval [12] should be extensible to multivariate quasi-ordinary 
polynomials. We give a new and more elementary proof of that fact and show 
how to compute parametrizations using our representation. We end with an open 
problem and an outlook in Section 5. In a short appendix we collect results from 
local commutative algebra for reference. 

Before we proceed we recall and fix some notions. Let E be a field of characteristic 
zero and X and Y integral E-schemes. By E(X) and E(F) we denote the respective 
function fields. A rational map tt : Y X is given by a tuple (V,t:) s.t. V <^ Y 
is open and tt : — > X is a regular morphism. Note that we do not restrict to 
schemes of finite type here. In particular all regular morphisms are rational maps. 
Two tuples (Vi, TTi) and (V2, 7^2) arc equivalent, or define the same rational map, if 

7ri|viny2 — t^^Wivw^.- 

Assume that two maps send the generic point oiY to j) & X (its image is always 
defined for rational maps). Then (14, tti) and (V2, "n-i) are equivalent iff the induced 
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inclusions of fields Ox,p/^x,p ^ E(y) are the same (where mx,p C Ox,p is the 
maximal ideal). In particular if tt is dense, i.e., p is the generic point of X, we get 
an inclusion E(X) ^ E{Y) determining tt. 

Note, however, that not all such field inclusions yield rational maps under our 
assumption since we have not yet restricted to schemes of finite type over E. E.g., 
let X := SpccE[a;], Y := SpecE[2;](^) and tt :Y ^ X he the morphism induced by 
localization. Then tt induces an isomorphism of function fields E(X) ~ K{Y). Nev- 
ertheless TT has no rational inverse. A rational map with inverse is called birational 
(or also a birational transformation). 

Further it is easy to see that dense rational maps may be composed. A rational 
map has a domain of definition, which is the maximal open set on which it can be 
defined (equivalently, the union of all such open sets). 

2. Formal Desingularizations 

For this section we denote by X and Y separated, integral schemes of finite type 
over E. Further we assume that both are of dimension n. All (rational) maps will 
be relative over E. 

Let {A,m) be a valuation ring of E(X) over E (where m is the maximal ideal). 
If A is discrete of rank 1 and the transcendence degree of A/m over E is n — 1 
then it is called a divisorial valuation ring of E{X) over E or a prime divisor of 
E(X) (see, e.g., [33, Dcf. 2.6]). It is an essentially finite, regular, local E-algebra of 
KruU-dimension 1 {i.e., the localization of a finitely generated E-algebra at a prime 
ideal, see [40, Thm. VI.14.31]). 

Let {A,m) be a divisorial valuation ring of E{X) over E. By [18, Lem. II. 4. 4.] 
the inclusion A C E(X) defines a unique morphism Spec QF{A) — > X and therefore 
a rational map Spec A X sending generic point to generic point. Composing 
this with the morphism obtained by the m-adic completion A ^ Awe get a rational 
map Spec A X. 

Definition 2.1 (Formal Prime Divisor). Let {A,m) be a divisorial valuation ring 
of E(X) over E. Assume that the rational map Spec A X (as above) is actually 
a morphism if : Spec A X {i.e., defined also at the closed point). Then ip is 
a representative for a class of schemes up to X-isomorphism. This class (and, by 
abuse of notation, any representative) will be called a formal prime divisor on X . 

Hence we may compose a representative ip with an isomorphism Spec B — > Spec A 
to get another representative for the same formal prime divisor. We have an iso- 
morphism A = WipftJ with := A/mA = A/m. In the sequel we will sometimes 
assume that A is already of this form, i.e., ip : SpecFj^p] — > X. The isomorphism 
is an instance of Cohen's Structure Theorem for regular rings, see Theorem A.l. 

Formal prime divisors provide an algorithmic way for dealing with divisorial 
valuations; A formal prime divisor yields an inclusion of function fields E(A') ^ 
¥ip {{t)). Vice versa, by what was said above, ip is determined by this inclusion. So 
it is this piece of information one has to represent. Composing this inclusion with 
the order function ordt : F^ {(t)) — » Z we get the corresponding divisorial valuation. 

We want to single out a special class of formal prime divisors. 

Definition 2.2 (Realized Formal Prime Divisors). Let p E X he a regular point 
of codimension 1. The formal prime divisor 

Spec O^p X 

(given by composing the canonic morphism Spec Ox,p X with the morphism 
induced by the completion Ox,p Ox,p) is called realized. 
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If X is normal then all generic points of closed subsets of codimension 1 are 
necessarily regular [18, Thm II.8.22A]. Therefore there is a one-one correspondence 
of realized formal prime divisors and prime Weil divisors. Another important fact 
is that we can compare formal prime divisors of birationally equivalent schemes 
under certain conditions. 

Lemma 2.3 (Fullback of Formal Frimc Divisors). Let X and Y he S-schemes with 
structure morphisms px '■ X S , py '■ Y ^ S and assume that py is proper. Let 
TT : Y X be a birational transformation of S-schemes (meaning that py = px'K 
as rational maps). 

Then each formal prime divisor : Spec A ^ X lifts uniquely to a formal prime 
divisor on Y, i.e., there is a unique formal prime divisor Tr*ip : Spec^ — s- Y s.t. 
tt(tt*ip) ~ if as rational maps. 



Proof. Consider the commuting diagram 
SpecQF(l) ^ SpccE(X) 



Spec A 



■SpecE(y) 




^Y 



where the vertical arrows are given by restricting identity maps to the germs at the 
generic points and the upper arrows are induced by ip and tt~^ respectively. The 
two squares trivially commute and the triangle commutes because tt was assumed 
to be a birational transformation of S'-schemes. 

Now assume that a lift Tr*(p as in the claim exists, then we must have 7r*(/3 = 7r~^(/? 
because tt is birational. This (a priori only rational) map would fit into the follow- 
ing contracted diagram: 



SpecQF(A) 



Y 



IT 



Spec A ■ 



S 



Now there exists a unique regular morphism ip : Spec A Y which fits into this dia- 
gram by applying the valuative criterion of properness to py (see [18, Thm. 114.7]). 
Since a rational map is uniquely determined by the inclusion of function fields we 
see that 7r*(f = tp. □ 

Corollary 2.4 (Fullback along Froper Morphisms). Let tt : Y X be a proper, 
birational morphism. A formal prime divisor on X lifts to a unique formal prime 
divisor on Y . Vice versa a formal prime divisor on Y extends to a unique formal 
prime divisor on X , hence tt* is a bisection. 



Proof. This is obtained by applying Lemma 2.3 to tt and tt ^ where S 
We will apply the operator tt* also to sets of formal prime divisors. 



X. 



□ 



Definition 2.5 (Center and Support). Let ip : Spec A X he a formal prime 
divisor. We define its center, in symbols center((/3), to be the image of the closed 
point. Further the support of a finite set of formal prime divisors S is defined as 
suppx('5) := {center(93) | ip G S}, i.e., the closure of the set of all centers. 

Example 2.6. Let X := Spec A with A := Q[x, y, z]/ {x"^ +y'^ — z^) be the cone over 
the circle and Y := SpccB with B := Q[x' , y', z']/{x'^ + y'^ - 1) the cylinder. The 
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strict transform tt under the blow up of the origin has an affine chart Y X given 
by the homomorphism A ^ B : x x'z' , y ^ y' z' , z i— > z' . The generic point of 
the exceptional divisor in Y is the prime ideal {z'). In this case we have a trivial 
isomorphism Bi^z') = QF(Q[.t', y']/ {x'-^ + y'^ — Then a formal prime divisor 

is induced by the homomorphism 

A ^ QF(Q[x',2/']/(^" + y" - : ^ ^ x'z',y^ y'z\z^ z' . 

We can compose this homomorphism with an arbitrary isomorphism of rings to get 
another representative of the same formal prime divisor. E.g.^ since the circle is a 
rational curve we can change the coefHcient field by QF(Q[x', y']/ {x'"^ + y'^ — 1)) 
Q(s) : x' I— > Y^^T, 2/' ^ "^"^.1 and map, say, z' ^ t + t^ + . . . . Now a formal prime 
divisor ip : SpecQ(s)[t] ^ X is induced by 

. 2s{t + t^ + ...) {-1 + S^){t + t^ + .. .) ^ o 

1 + 1 + 

One finds ccnterx('y3) ~ {x, y, z) (which is the preimage of the prime ideal (t)) and 
by construction we know that TT*(p is realized with center (z') in Y . 

Now we are in the situation to define formal dcsingularizations. 

Definition 2.7 (Formal Description of a Desingularization). Let tt : y ^ X be a 

desingularization, i.e., tt is proper, birational and Y is regular. Let 5 be a finite 
set of formal prime divisors on X . We say that 5 is a formal description of tt iff 

(1) all divisors in ■n*S are realized, 

(2) 7r~^(suppjf (iS)) = suppy(7r*5) and 

(3) the induced morphism y \suppy(7r*iS) X \ suppj^j- (5) is an isomorphism. 

The set S itself consists of formal prime divisors on X and makes no reference 
to the morphism tt. By another definition we can avoid mentioning an explicit tt. 

Definition 2.8 (Formal Desingularization). Let iS be a finite set of formal prime 
divisors on X. Then S is called a formal desingularization of X iff there exists 
some desingularization tt s.t. 5 is a formal description of it. 

Informally speaking the set S makes it possible to treat divisors on Y effectively, 
although we haven't explicitly represented y as a whole; Indeed, realized formal 
prime divisors correspond bijectively to usual prime divisors on the regular scheme 
y. The set of formal prime divisors now divides into two classes: those within 
■K*S and those with center in y \ suppy (7r*iS). The latter can be dealt with on the 
isomorphic scheme X \ supply (iS). Therefore formal descriptions are an appropriate 
algorithmic tool to work with invertible sheaves on Y . 

In the case of surfaces it is easy to see that the existentially quantified tt in the 
above definition is actually unique up to isomorphism. Therefore S really identifies 
a desingularization in the common sense. Vice versa, every desingularization can 
be described formally by completing the stalks along the exceptional divisors. 

Theorem 2.9 (Uniqueness of Surface Desingularization). Let S he a formal desin- 
gularization of X . If -Ki : Yi ^ X and ^2 ■ Y2 ~* X are two dcsingularizations 
described by S (in the sense of Definition 2.7) then Y\ and Y2 are isomorphic as 
X-schemes. 

Proof. Let 6'"^ : — > X be a minimal desingularization (see [26, Cor. 27.3]) of 
X and let 5*^°-' be obtained from (&(°))*cS by subtracting all realized formal prime 
divisors. Then tti and tt2 factor through X^^') yielding a commuting diagram 
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and ttJ" and f-re both described by S^^^ . 

Moreover these maps are proper, birational morphisms between regular surfaces, 
so they factor into a finite sequence of point blow ups (see [18, Cor. V.5.4] which 
holds also for the case of a non-closed ground field). The set of possible centers 
for the first blow up is exactly supp;5(-(o) Choose a center, compute the blow 

up 6(1) : and set S^'^^ to be the set excluding the unique 

formal prime divisor that is turned realized and centered along the exceptional 
divisor. Again tt^^ and tt^^ factor through morphisms ttJ^'' : Yi — > X^^^ and 
7i'2^^ : Y2 —> X^-^) described by S'^^K Going on hke this and setting I := #<S(°\ we 
find that ttJ''' and ttj''' are described by 5'-'' = and hence are isomorphisms by 
Definition 2.7. □ 

It is not so clear whether a similar statement holds in higher dimensions when 
minimal resolutions are not available. 

Remark 2.10 (Formal Desingularization of Reduced Schemes). The above defini- 
tions can and will be used in a more general setting. Namely, if X = IJ^ is 
the decomposition of a reduced, equidimensional (no longer integral) scheme into 
irreducible closed subschemes then any morphism SpecF|i| ^ X is actually a mor- 
phism to one of the Xi. We call it a formal prime divisor if the corresponding 
SpecF|<] Xi is a formal prime divisor. We call it again realized iff it corre- 
sponds to the completion at the germ of a regular codimension 1 point in X (not in 
Xi\). Lemma 2.3 and Corollary 2.4 remain valid in this setting (where a birational 
morphism between two reduced schemes is a morphism that induces birational 
morphisms on all irreducible components). The definitions of center, support and 
formal desingularization carry over straight forward. 

3. The Method of Jung Revisited 

In this section we describe the desingularization of surfaces after Jung. 

3.1. Theory of the Method. First we will view a projective surface as a certain 
covering of a smooth surface. Then we modify the covering such that after passing to 
the integral closure the remaining singularities are very simple and can be resolved 
by point blow ups. 

3.1.1. Projective Surfaces as Ramified Coverings of the Plane. Consider a projec- 
tive hypersurface X C P|. We want to view X in a slightly different way. Let 
Po G P| \ X be a closed rational point (i.e., its residue field is isomorphic to E) 
and TT : W V the linear projection from po where W Pg \ po and V := Pg. 
This projection defines a line bundle. Its restriction tt jx^ — > ^ is a Noether 
normalization of X, i.e., a finite morphism onto the projective plane. We subsume 
the governing properties in the following definition: 

Definition 3.1 (Ramified Coverings). Let tt : W ^ V he a line bundle s.t. F is a 
regular, integral surface over E. Further let X C Vl^ be a reduced hypersurface s.t. 
tt\x '■ X V is finite. The tuple {tt,X) is then called a ramified covering. 

That is our notion of ramified covering comprises that we are dealing with sur- 
faces and that the covering surface is embedded in a line bundle over the base. Now 



FORMAL DESINGULARIZATION OF SURFACES 



7 




Figure 1 . Embedded Desingularization of the Discriminant 

The left side shows the local picture of the example surface X together with its discriminant curve. On the right 
we depict a chart Xq obtained from the embedded desingularization of the discriminant curve. 

we want to define the ramification locus of such a covering, i.e., the locus where 
the covering 7r|jf is "not locally trivial" (more precisely, not etale). Since tt is a line 
bundle and Tr\x is finite, we can find a covering {Ui}i of V by affine open subsets 
s.t. 0\Y{'^~'^{Ui)) = Ov{Ui)[z] is a polynomial ring in one variable and X is given 
by a monic, squarefree polynomial fi G Ov{Ui)[z]. 

Definition 3.2 (Discriminant Curve). Let {n : W V,X) be a ramified covering. 
Let {Ui}i be a covering of V as above and fi G Ov{Ui)[z] a set of polynomials 
defining X. The ideals (disCz(/i)} C Ov{Ui) define an invertible sheaf of ideals and 
the corresponding subscheme is called the discriminant curve -D,r|x — ^■ 

The discriminant curve is actually a concept independent of the concrete covering 
and embedding of X into W. It depends only on tt\x- Further the covering is locally 
trivial except over i?7r|x- 

Example 3.3. Consider the surface X C Fq given by F = with F := Xq + 
3xqX2X3 + Xqx\x2 + ?>x'^x^x\ + X2x\ £ 'Q[xo, Xi,X2,x^]. Sincc F is monic in xq the 
surface X doesn't contain the point p := (1 : : : 0). The line bundle defined 
by the projection of X \ {p} to the plane = is given by dchomogcnizing with 
respect to xi, X2 and x^. 

In the last chart, the defining equation has the form / := + 31'^"* + v?vur' + 
"iv^vo^ -\- v'^ £ 'Q[u,v\[w\ where wc have mapped xq ^ w, xi ^ u, X2 ^ v and 
X3 1-^ 1. We have the local discriminant disCu,(/) = 12Qu^v^^{u'^ — 64ti). Figure 1 
left displays the surface. 

3.1.2. Embedded Desingularization of the Discriminant Curve. We would like to 
give a more detailed study of the covering, but for general ramified coverings (tt, X) 
this is hard. The complexity of X resp. of the covering map is partially reflected in 
the discriminant curve ^'ttIx ■ Recall that a closed point p is called normal crossing 
for the embedded curve D^\^ C y if the curve is locally defined by d^.^d'y' u S Ov,p 
where u is a unit, {dx,dy} is a local system of parameters and ex,ey > 0. The 
whole curve is called normal crossing if it is normal crossing at every closed point. 

Definition 3.4 (Nicely Ramified Coverings). Let (tTjX) be a ramified covering. If 
Z?7r|j^ is normal crossing then we call (tt, X) nicely ramified. 

We can always modify a ramified covering to become nicely ramified. If we 
speak of normal crossing singularities of the discriminant curve of a nicely ramified 
covering, we mean the closed points where two components intersect. 
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Lemma 3.5 (Simplification of Coverings). Let {tt : W ^ V,X) be a ramified 
covering. There is a proper, birational morphism p : V' ^ V s.t. the ramified 
covering (tt' : W V',X') is nicely ramified (where W := W xy V and X' := 
X Xw W ). Further p is given by a finite succession of blow ups in closed points. 

Proof. The theorem on embedded desingularization of curves (see [24, Thm. 1.47]) 
shows the existence of the morphism p; Indeed it says that after a finite number 
of blow ups in the singular points of the reduced curve the puUback of Djrlx 
normal crossing. For showing the lemma it remains to prove that constructing the 
discriminant curve commutes with base extension. This is left to the reader. □ 

Example 3.3 (continued). The curve defined by disc^(/) has a complicated singu- 
larity at the origin which needs to be resolved. One of the chart maps is given by 
u I— > uv, V 1-^ u^v^ and transforms disCu,(/) to 729m'^^?;*^(m^w — 64) which describes 
a curve with a normal crossing intersection at the origin. 

The embedded desingularization of the discriminant curve can be applied to the 
surface X (by mapping w i— > w) to obtain a surface with chart Xq. We obtain the 
new local equation fa ^ + Zu^v^w^ + u^v^w^ + Zu'^v^vP + u^v^ with disc^j (/o) = 
729u'^^t;'*''(w^u — 64), see Figure 1 right. 

3.1.3. Desingularization of Toroidal Surface Singularities. The structure of nicely 
ramified coverings depends crucially on the celebrated Theorem of Jung-Abhyankar 
(see Theorem A. 5). A polynomial fulfilling the conditions of the theorem is called 
quasi- ordinary. In its original form the theorem is not precise enough for our 
purposes, for example, because the statement doesn't involve the coefficient fields 
of the power series solutions. In the sequel IE|x'"] will denote the ring of power 
series with coefficients in a field E, variables xi^ . . . ,Xn and exponents in F H M"q 
(the non-negative orthant of a full rational lattice) . We need the concept of rational 
parametrizations (for a refined version see Definition 4.9): 

Definition 3.6 (Parametrizations). Let / G E|a;i , . . . , x„| [z] be quasi-ordinary. 
I.e., f is monic, squarefree and s.t. disCz{f) ~ x']^ ■ ■ ■ x'^u{x) where u(0, . . . , 0) ^ 0. 

We call (cr, a) with a S Aut(E'|a;i, . . . , Xn\ \ E') and a S E'lx'"! a parametriza- 
tion of / if E C E', Z" C F and crT^(/)(a) = 0. Let g\f be an irreducible 
factor s.t. (j^^{g){a) = 0. We call (cr, a) rational if the induced homomorphism 
IC(E|.Ti, . . . ,x„|[z]/(5)) E'fa;'"] which maps z i— > a and 7 1-^ a{j) for 7 e 
E|a;i, . . . ,Xn\ is an isomorphism. A set of rational parametrizations for / is called 
complete if it is in bijective correspondence with the irreducible factors of /. 

Here we used a^^ for the lifting of the automorphism to the polynomial ring by 
coefficient-wise application. We will show later, that we can actually compute such 
rational parametrizations. 

Theorem 3.7 (Existence of Rational Parametrizations). Let f e E|a;i, . . . ,a;„][z] 
be a quasi-ordinary polynomial. Then a complete set of rational parametrizations 
of f exists and can be computed. Moreover, if actually disc2(/) = x^^ ■ ■ -x^u^x) 
where m < n and w(0, ... ,0) 7^ then the exponent lattices of the power series 
rings will be of the form F x Z"~™ for some m-dimensional rational lattice F. 

Proof. For the first statement see Algorithm 11 and Corollary 4.15 of Section 4.1. 
In Lemma 4.14 the relation between the computed parametrizations and the frac- 
tionary power series roots of / is explored. Together with Theorem A. 5 this gives 
the statement about the exponent lattice. □ 

Example 3.3 (continued). The transformed polynomial /o G Q[m, is quasi- 
ordinary. In this case a complete set of rational parametrizations is given by only 
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a single parametrization (cr, a) with cr : Q|u, u] — >■ Q|u, u] : u > — 8u, w i— > — w and 

69 810 1011 1413-|1815 ,^2217 -,2619 

a := —8u6 V 6 +8u6 V s —Au^vS +u6ti6 -~^u e v 6 +^it6u6 ^-^u 6 t) 6 +. . . 

Wc will sec later that the simple form of a is no coincidence. We further have 
a e with T := Z(0, i) + Z(i, i). This lattice is shown in Figure 2 below. 

Complete sets of rational parametrizations describe very explicitely the structure 
of the integral closure of E|a;i, . . . , a;„][z]/ (/). 

Lemma 3.8 (Decomposition Induced by Rational Parametrizations). Let f S 

El^i, . . . ,a;„|[z] be a quasi-ordinary polynomial andV = {(cri,Q!i)}i with {ai,ai) S 
Aut(Ei|a;i, . . . ,Xn\ \ E;) x E^Js.'"'] a complete set of rational parametrizations of f. 
Then the homomorphisms 

ipi : Elxi, . . . ,Xnl[z]/{f) Eil^^'l : 7 ^ '^ii.l) for 7 G E|[a;i, . . . ,a;„] and z a,; 
can be composed to a homomorphism 

i> : E[.Ti, . . . ,x„l[z]/(/) ^ n^'fe^'l : a ^ (. . . , ■ • ■ ) 

i 

which lifts to an isomorphism -0 : IC(E|xi, . . . , Ili ^ife'"'!- 

Sketch of Proof. By the definition of rational parametrizations it is enough to show 

IC(EIxi, . . . , XnM/if)) - n. IC(E[a:i, . . . , Xn\[z]/{f^)) 

where the righthand side runs over all irreducible factors fi of /. By [13, Exer. 
2.26] this can be shown using orthogonal idempotents. More precisely, if / = /11/12 
where hi and ft.2 are factors without common divisor then set ei := /i2/(^i + ^2) and 
62 := hi/{hi + h2) to be elements in QF(E|2;i, . . . , a;„][z]/ (/)). (For this one checks 
that the denominators are no zero-divisors.) Then one easily computes ei 4- 62 = 1, 
6162 = and Ci = ei{ei + e-i) = ef + 6162 = e|. The idempotency relations imply 
in particular that 61,62 G IC(E|xi, . . . , a;„| [z] /(/)). So the integral closure splits: 

IC(E[a:i, . . . ,x„I[z]/(/)) - 61 IC(EIxi, . . . , a;„l[z]/(/))x62lC(E[a;i, . . . ,xJ[z]/(/)) 

Finally one shows 6jC(Ela;i, . . . , a;„l[z]/(/)) = IC(E[a;i, . . .,XnM/{h^))■ □ 

Corollary 3.9 (Singularities of Ramified Coverings). Let {7t,X) be a ramified cov- 
ering, V : X X the normalization morphism. Then the (isolated) singular points 
of X lie over the singularities of the (reduced) discriminant curve. Moreover, over 
normal crossing singularities the singularities of X are toroidal. 

Proof. Let g G X be a closed point and p := TT{v{q)). Assume that p either does 
not lie at all on -DttIa-i ^i'-'^ '^'i ^ regular point of the reduced curve, or is a 
normal crossing singularity. In each of these cases Dt^\^ can locally be defined by 
d^'^dy^u G Ov,p where {dj;,dy} is a local system of parameters and e^; > 0, 6^ > 0. 
There is an isomorphism between the completion of Ov,p and F|a;,y] (where F is 
the residue field of Ov,p) which maps dx >^ x and dy 1-^ y, see Theorem A.l. We 
can as well assume that Ov,p =^lx,y}i, d^ = x and dy = y. 

The completion of the fiber of 7r|x can be defined by the vanishing of a polynomial 
/ G FJx, which is quasi-ordinary. Building the integral closure commutes with 
completion, see Lemma A. 3. Then Lemma 3.8 implies that the completion of ^ 
is isomorphic to a power series ring with fractionary exponents. But those are the 
completions of the distinguished stalks of affine toric surfaces, hence, all such points 
q can at most be toroidal. 

If p does not lie on i?,r|x then e^. = Cy = 0. If it is a regular point of I?7r|x then, 
say, 62; > and Cy ~ 0. By Theorem 3.7, in both cases, the power series rings in 
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Figure 2. Lattice and Dual Lattice 

We show the non-negative quadrant of a lattice F and its dual (i.e., the set of all linear forms to Z identified 
via scalar multiplication of vectors). In the dual we have inscribed an minimal sequence of generators ni, . . . , n^. 
Indices are chosen s.t. subsequent pairs correspond to neighboring vectors. 

the parametrizations have exponent lattices of the form (-Z) x Z. These rings are 
regular and by faithful flatness ^ must be regular itself. □ 

Example 3.3 (continued). For our example this means that the integral closure of 
Qfu, v}[w]/ (/o) is isomorphic to Q|(u, v)^}. This fits well to the picture in Figure 1 
that suggests that the surface Xq, though singular, has only a single analytic branch 
at the origin. 

Modifying a ramified covering to become nicely ramified is constructive, since 
embedded desingularization of curves is. Passing to the normalization of a scheme 
of finite type over E is constructive, since computing the integral closure of finitely 
generated E-algebras is (see, e.g., [11]). By the above corollary we are left with the 
task of desingularizing normal toroidal surface singularities. 

These singularities have first been studied by Jung [22]. As noted in [9, Lcct. 2, 
§2] a normal toric surface can be desingularized by a finite sequence of point blow 
ups. This property can be transferred to toroidal singularities by Lemma A. 4. So 
we could obtain a desingularization by computing a sequence of point blow ups of 
the normalization X and, hence, Jung's method is already constructive. 

However, computing normalizations is not cheap. Lemma 3.8 shows that rational 
parametrizations anyway describe the integral closure. In the following section we 
will therefore follow a different approach, always with formal desingularizations in 
mind. We will also benefit in other places from the additional flexibility obtained 
by applying formal isomorphisms. 

The final lemma of this section has its origins in Hirzebruch [21] (see also [25]) 
who first gave explicit formulas for the desingularization of toric surface singularities 
using continued fractions. Afterwards the arrival of toric geometry (c/. [15]) has 
introduced new terminology and means of description. Recall that the dual of a 
fuU lattice F C Q" is the set of all linear forms n G Q" with n(m) e Z for all m £ F. 

Lemma 3.10 (Formal Desingularization of Toric Surfaces). Let SpecE[^^] for 
some rational lattice F C be a toric surface and ni,...,n; G F^ fl R>q an 
ordered and minimal sequence of monoid generators as in Figure 2. Then the set of 
morphisms tpi : SpecE(s)[[i| SpecE[2.'"] given by the E-algebra homomorphisms 

E[x^] E(s)p]] : ^ 5n,(m)^n.+i(m) 

for l<i<l^2isa formal desingularization of SpccE[x^]. 

Proof. Desingularizations of toric schemes can be constructed using special fans 
in the dual cone and lattice F^ n M>q (see, e.g., [15, 10]). Applying this con- 
struction to the fan whose one-dimensional faces are given by the rii one obtains 
a desingularization tt. In fact the morphisms tt^ : SpecE[s,t] SpecE[x^] given 
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by E[^^] E[s,t] : x"' ^ gn,{m)^n,+i{m) for 1 < i < Z - 1 are isomorphic to the 
restriction of tt to an open covering by affine charts. 

The exceptional locus {i.e., the 7r-preimage of the isolated singularity) is a finite 
union of divisors. We get a formal description of tt by completing the local rings 
along these divisors. They are given, for example, by the prime ideals (t) C E[.s, t] 
in each of the charts 1, . . . , Z — 2. For the completions of the local rings we then 
have trivially IE[s^i) ^ IE(s)[[t]. □ 

Remark 3.11 (Minimal Desingularization by Point Blow Ups). The desingulariza- 
tion described by the above set of morphisms is the minimal one obtained (up to 
isomorphism) by subsequently blowing up isolated singular points. To show this 
one observes that every single such blow up gives (non-afiine) toric surfaces that 
are described by fans in n R>q. Elementary arguments about the lattice and 
its dual show that the one-dimensional faces of these lattices are always given by 
one of the vectors . On the other hand any fan defined by a proper subset of the 
generators is associated to a still singular surface. 

Example 3.3 (continued). This lemma can be applied as follows. The normal- 
ization Xq of Xq has an isolated singularity which is formally isomorphic to the 
distinguished germ of a toric surface with coordinate ring Q[(u,w)'"]. 

The lemma says that this toric surface is formally desingularized by mapping 
(m,u)" ^ sni(m)^n.+i(m) ^ Q(s)|t] for i £ {1,2} with as in Figure 2. Let's look 
at this for i = 1 in terms of algebra generators v^^^ jU^^'^v^^^ ,u S Q[(m, w)^]. We 
have to map 

„1 ^ g((0,6),(0,l/2)>^((l,4), (0,1/2)) ^ ^3^2^ 

Jy^ ^ ^((0.6),(l/3,l/6))i((l,4),(l/3,l/6)) ^ 
^^S<(0,6),(1,0))^((1.4),(1,0)) 

Composing with the homomorphism induced by the rational parametrization 
{a, a) we get homomorphisms ipi : Q[u,v][w]/ (fo) — *■ Q(s)[t]: 

^,:u^ -8s% V ^ -t\ w ^ -Ss'^f + 8/t« - is^°f + s^H^' - ^s^^'t^^ + ..., 

(P2 : u ^ -8s^^ V ^ -sH^, w ^ -8sV^ + Bs^t^ - is^f + s"<^ - ^s^^^" + ... 

The two induced morphisms SpecQ(s)|i| —>■ Xq are the formal prime divisors 
centered at the origin that become realized on a desingularization of Xq by point 
blow ups. 

But this is not yet a formal desingularization of Xq. As can be seen in Figure 1, 
the singular locus of Xq is the union of the u-axis and the ii-axis. So there should 
also be formal prime divisors supported on these lines. Let's first consider u = 0. 
By mapping t and w s we transform /o to + ?>s^t^w'^ + s^t^w^ + is^f^w^ + 
s^t^ G Q(s)[t][u']. Now we compute a univariate rational parametrization (cr, a): 

a:t^ ^t,a^ -^^th - ^tl + ^^tl + -^^tl + ^^t^ + • • • G Q{s) [7] [t^i/^)^] 

Here 7 is an algebraic element satisfying 7^ + = 0. Canceling exponent denomi- 
nators this yields a homomorphism (fQ : Q[u,v][w]/ (/q) Q(s)[7]|t]: 



ipQ : ui 



This homomorphism corresponds to completing the germ at a generic point in Xq 
above u = by Lemma A. 3. This germ is not modified by subsequent point blow 
ups. So we get a further formal prime divisor. A last one is supported on v = 0: 

ip3:u^s,v^ -64s^°^^ w h-> -512s^^i^ + 512s^H^" - 2565^"^" + 645^"^*^^ + . . . 
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Note that this procedure exphcitly produces the residue fields at the generic 
points of the exceptional divisors in a desingularization. They correspond to the 
coefficient fields of the power series. Finally, composing all four formal prime divi- 
sors with the morphism Xq — > X which was obtained by the desingularization of 
the discriminant curve, we get a part of a formal desingularization of X. 

3.2. A Divide and Conquer Approach. Now we want to cast the theory of the 
previous paragraph into explicit algorithms. We want to give a formal description 
of a Jung desingularization of a ramified covering {n : W ^ V, X). Note that such a 
desingularization is always relative to an embedded desingularization p : Vp ^ V of 
the discriminant curve D^^ ,^ . Because then {iTp : Wp Vp, Xp) for Wp := W XyVp 
and Xp := X X\Y Wp is a nicely ramified covering and we can define the Jung 
desingularization Up : Yp X to he obtained by normalizing Xp and successively 
blowing up singular points, see Remark 3.11. 

If we wanted to avoid redundant blow ups, we could already fix the embed- 
ded desingularization p to be minimal, i.e., obtained by blowing up a point of the 
discriminant curve only when it is a non-normal crossing singularity. For compu- 
tational reasons, we choose a slightly different desingularization, see Remark 3.17 
below. Also, the definition of a formal description S of Hp leaves a bit of choice for 
S, e.g., one may always add formal prime divisors which are realized on the desin- 
gularization. We get rid of this ambiguity by requiring that supp y (5) = tt" ^ {D^ | ^ ) 
and say we compute a formal description of Hp above i'lrix- 

Remark 3.12 (Divide and Conquer Paradigm). We will dcsingularizc the discrim- 
inant curve and compute the formal prime divisors for the surface at the same 
time. The substitutions involved in computing the curve desingularization make 
the defining equations more complicated. Therefore our paradigm must be to com- 
pute formal prime divisors "as early" as possible. In other words, if we know that 
the surface (resp. its normalization) remains unchanged in a subset (up to isomor- 
phism) by further blow ups of the discriminant curve we immediately compute the 
formal prime divisors centered in that set. 

Let p : Vp V he an embedded desingularization of the discriminant curve. 
Our divide and conquer approach (in particular Algorithm 3 below) relies on the 
following facts: 

• Let p G i?7r|x be a point s.t. p is not an isomorphism at p. Let po ' 
V ^ V he the blow up at p and (tt' : W' — > V',X') the induced ramified 
covering. Then p factors as p = piPo where pi '. Vp V' is an embedded 
desingularization of D^,^^,, Xp X factors through X' X and also 
7r'~-'^(D^/|^, ) must be equal to the support of the puUback of S. So it is 
equivalent to compute the formal description of IIp^ above D^^i^^,. 

• Now Po is an isomorphism everywhere except at p. We can split the compu- 
tation of the formal desingularization Tip into two parts; the computation 
of formal prime divisors which are not centered above p on the one hand 
and those which are centered above p, or equivalently, whose puUbacks (see 
Corollary 2.4) are centered above the exceptional divisor p^ ^(p). Comput- 
ing the latter will be delegated to a recursive call. 

• When blowing up the (not necessarily rational) point p we may first apply 
a morphism to the projection plane that induces a formal isomorphism at 
p because of Lemma A. 3 and Lemma A. 4, compare Remark 3.14 below. 

Finally, Example 3.3 has shown that in the case of nicely ramified coverings 
we have to compute formal prime divisors in two ways: Those which are centered 
above the components of the discriminant curve are obtained using Lemma 3.8 with 
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n = 1 , and those which are centered above the normal crossings of the discriminant 
cmve using a combination of Lemma 3.8 with n = 2, Lemma 3.10 and Lemma A. 4. 

3.3. The Algorithm. In the foUowing algorithmic descriptions we allow subsets 
of a set A, which in our notation will be elements of 2'^, as data types. These 
will either be finite sets or they will be finitely generated ideals of a ring A. So 
it is clear that they can be represented. For simplicity of notation, we also allow 
passing of homomorphisms from polynomial rings in a finite number of variables. 
They can obviously be represented by the images of their generators. We assume 
that we can represent power scries which will be explained later in Section 4.1. If 
: A — > i? is a homomorphism of rings, we write again (p^^ : A[w\ — > B[w\ for 
its lifting to the corresponding polynomial rings by coefhcient-wise application and 
map^ : 2"* 2^ for the function on subsets defined by element-wise application. 

Let X C P| be a closed hypersurface. Following Section 3.1.1, we first have to 
produce a ramified covering. This is done in algorithm DesingGlobal. By E (resp. 
E(s)) we denote the algebraic closure of E (resp. of the rational function field). 
If it shows up in the return type of a signature we actually mean that the result 
involves some finite field extension (of transcendence degree 1), i.e., we do not rely 
on a system for computing with algebraic closures. 

The algorithm will be formulated for reduced surfaces, see Remark 2.10. This 
is due to the fact that we will transform our surface by extending the base field in 
certain steps and cannot assure that the transformed surface remains integral, even 
when the original surface was, see Remark 3.14 below. 



Algorithm 1 DesingGlobal(F : E[xo, . . . ,2:3]) : 2E[^o,...,^3]-e(s) W 
Require: A squarefree homogeneous polynomial F ^ 0. 

Ensure: A finite set of homomorphisms E[a;o, . . . , X3] Fp| factor- 

ing through E[xo, ■ ■ ■ ,X3]/ (F) s.t. the induced morphisms SpecF|t| 
Proj E[a;o, . . . , X3]/ (F) are a formal desingularization. 

1: Let (j) : E[a;o, . . . , xa] E[.to, . . . , xa] be a linear automorphism s.t. 

0(F)(1, 0,0, 0)7^0; 

2: ipi : E[xo, . . . , xa] — > E[m, v][w] : xq i— > w, xi i~> 1, X2 ^ u, 2:3 t-^ v; 

3: S := map^^^^^^DesingLocal(V'i0(F), (0)); 

4: V2 ■ E[2^0j ■ ■ • , X3] — > E[w, v][w] Xq 1-^ Xi 1-^ V, X2 i— > 1, ^3 i— > U] 

5: S := S\J mapp^p^^^DesingLocal(-020(-F), {v))] 

6: V3 ■ E[xo, . . . , x^] E[m, v\[w] : xq ^ w,xi ^ u, X2 I— > w, 2:3 i~> 1; 

7: S := SU mapp^p^,3^DesingLocal(i/'3(/)(F), 

8: return S; 



In line 1 we choose a linear automorphism of P| (represented by </)) s.t. the 
preimage of X under this automorphism is Noether normalized by a projection 
onto the plane xq ~ 0- 

Remark 3.13. The automorphism maps (1,0,0,0) to a point p ^ X. In an actual 
implementation one should find p s.t. most of its coordinates are zero and the rest 
are small integers. This preserves sparsity in 4>{F) and keeps coefficients small. 

For convenience of description, we will actually assume that (j) is the identity. 
Then the ramified covering (tt : W — > V,X) is given as follows: We set W := 
'^E \ {(1 : : : 0)}, V is the plane xq = 0, tt is the corresponding linear projection 
and X is defined by the vanishing of F monic in xq. 

The algorithm produces a set of homomorphisms S representing a formal desin- 
gularization. Therefore we cover V by open subsets Xi ^ (given by the ■0,) for 
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1 < i < 3 and call algorithm DesingLocal for each of those in lines 2 to 7. The 
latter algorithm produces formal desingularizations of the respective affine subsets. 
Because of the huge overlaps we add focus ideals to each call. 



Algorithm 2 DesingLocal(/ : E[u,v][wlT : 2^^"'"^) : 2E[",'']M^E(s)ltl 

Require: A squarefrec polynomial f ^ 0, monic in w, and a focus ideal JF. 
Ensure: A finite set of homomorphisms E[u, F|t] factoring through 

E[u, w] s.t. the induced morphisms SpccF[[t| SpccE[w, i'][w]/(/) form 

the subset of a formal desingularization which is centered above the closed 
subset defined by T. 
1: d := SquareFreePart(disCtu(/)); E :~ IrredFactors(c?); 

E with J^cie) DivisorsAboveCurve(/, e); 
3: for {uo,vo) e ZeroSet(jr+ {d,dd/du,dd/dv)) do 
4: -0 : E[m, v] ^ E'[u, w] : u i-^ u + uq, V V + vq; 

5: S := SU mapp^p^T»DesiiigRec(V;T"'(/), {^/.(e) | e G £; and 7/'(e)(0, 0) = 0}); 
6: return S; 



In line 1 we compute the defining equation d of the reduced discriminant curve 
and its factors E. For all e S ii^ the prime ideal (e) corresponds to the generic point 
p £ of an irreducible component of the discriminant curve. 

No matter how the desingularization p : Vp V oi the discriminant curve looks 
like, it is always a succession of point blow ups. Therefore p will be isomorphically 
contained in Vp. The same holds for points in the normalization of X lying above p. 
According to our paradigm in Remark 3.12 we compute the formal prime divisors 
centered above these p already at this stage of the algorithm, see line 2, by calling 
algorithm DivisorsAboveCurve (see Algorithm 4) for each e G E which is in focus. 

In line 3 we compute the singular locus of the reduced discriminant curve, more 
precisely, that part which is in focus. For each of its closed points p we want to find 
the formal prime divisors centered above p. Therefore we apply the homomorphism 
tp of line 4 (which corresponds to moving p to the origin) and then call algorithm 
DesingRecursive in line 5. 

Remark 3.14 (Exploiting Formal Isomorphisms). At this point we would like to 
mention that ip may involve an algebraic field extension, namely, by the residue 
field E' of p. Therefore tp, which looks like a mere translation, is not an isomor- 
phism. Nevertheless the induced morphism of schemes SpecE'[M,T;] — > SpecE[M, ii] 
is formally isomorphic at p by Corollary A. 2. (Note that we can choose a system 
/i, /2 e E[m, v] of generators of the maximal ideal corresponding to p s.t. /i e E[u] 
is irreducible and the image of /2 in (E[u]/(/i ))[?;] is also irreducible.) But for 
computing a formal desingularization, we may well pass to a formally isomorphic 
scheme as a consequence of Lemma A. 3 and Lemma A. 4. Note that the introduc- 
tion of a field extension may further split the defining equation, i.e., ip^'^{f) may 
be reducible even though / is not. Also the discriminant factors might split again. 
This is another (more important) reason for computing the formal prime divisors 
above the components of the discriminant already in line 2. 

Remark 3.15 (Auxiliary Functions). The above algorithm depends on a couple of 
auxiliary functions which we are not giving in detail, their names are mainly self- 
explanatory: SquareFreePart should compute the squarefree part of a polynomial 
and IrredFactors is supposed to produce the set of irreducible factors of a poly- 
nomial. A comment on ZeroSet is in order. It expects a zero-dimensional ideal 
T C E[a;i, . . . ,Xn]- It should return a finite set of rt-tuples s.t. for each maximal 
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ideal f containing T there is exactly one tuple (^i, . . . G (E')" s.t. E C E' and 
the induced homomorphism E[a;i, . . . , a;„] — > E' : Xi i-^ lifts to an isomorphism 
from E[xi, . . . ,x„]/f. 

Remark 3.16. Computing the squarefree part and factorization was done in Al- 
gorithm 2 to keep the number of parameters small. Of course it would fit more 
naturally in Algorithm 1 to avoid multiple computations. 

The recursive algorithm DesingRecursive now implements the divide and con- 
quer paradigm, see again Remark 3.12. 



Algorithm 3 DesingRecursive(/ ; E[u,v][w],E : 2'^["-^l) : 2E[«.''lM^E(s)ltl 

Require: A squarefree polynomial f ^ 0, monic in w, and a set of polynomials E 
s.t. e(0,0) = for all e G i? and JleGB ^ ^^'^ squarefree part of disct„(/) in 
the local ring E,[u,v](^u,v)- 
Ensure: A finite set of homomorphisms E[u, F|t] factoring through 

E[u,?;][w]/(/) s.t. the induced morphisms SpecF|i| — > SpecE[M, form 
the subset of a formal desingularization which is centered above the origin. 
1: if IsNormalCrossing(£') then 
2: return DivisorsAboveCrossing(/, £■); 

3: (t>u ■ v] — > E[u, v] : w i—> uv, v ^-^ v; (pr, '■ E[u, v\ E[u, v\ : u i—f v,v i—>- uv] 
4: /„ := ^J^^if); fv := ^Jn/); 

5: Eu ■= mapg^^^(-g)/^ord(e)(£'); E^ := map^^^^ (gj/^ordce) (£'); 

6: S mapp^p0^T™DivisorsAboveCurve(/u, f ); 

7: for (uo,Wo) G ZeroSet((neg£;„ e,t;)) do 

8: ip : E[m, v] — > E'[it, w] : It > It + mq, v t-^ v + vq; 

9: 5 := 5 U mapp^p,0T™0„T™DesingRecursive('0T'"(/^)^ 

{?/'(e) \ eeEu and V'(e)(0, 0) = 0} U {v}); 
10: if exists e € E^ s.t. e(0,0) = then 
11: S := 5 U mapp^p0^Tu,DesingRecursive(/t,, 

{e\e& Ey and e(0, 0) = 0} U {v}); 

12: return S"; 



In line 1 we check, whether the origin of the reduced discriminant curve is a 
normal crossing. If this is the case we know (see Corollary 3.9) that the points 
above the origin in the normalization are toroidal singularities that can be desingu- 
larized by a succession of blow ups in isolated singular points. The corresponding 
formal prime divisors are computed by algorithm DivisorsAboveCrossing (see 
Algorithm 6) and returned. 

Otherwise we have to modify the discriminant curve by blowing up the origin. 
The two charts of the blow up are computed in lines 3 to 5: We determine the defin- 
ing equations /„ and /„ of the transformed surface and also modify the discriminant 
factors accordingly. 

Note that the homomorphisms are such that the exceptional divisor is given by 
V = in both charts. Over this exceptional divisor there are centered a couple of 
formal prime divisors. By the same reasoning as for Algorithm 2 (see Remark 3.14) 
we compute them immediately by calling DivisorsAboveCurve in line 6. 

Now we have to consider the points on the exceptional divisor which are singular. 
In one of the charts they are given by the set of line 7. As in Algorithm 2 we move 
these points to the origin and go into recursion, see lines 8 and 9. Now there is 
possibly one singular point left to consider, namely, the origin of the other chart. 
It is treated in line 11. 
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To complete the algorithm, it remains to show how to compute the formal prime 
divisors which are centered above irreducible components or normal crossings of 
the reduced discriminant curve. 

3.3.1. Divisors above Generic Points of the Discriminant Curve. We now give an 
algorithm for computing the formal prime divisors which are centered above generic 
points of the discriminant curve. This is easy by what we have developed so far. 



Algorithm 4 DivisorsAboveCurve(/ : E[u,v][w], e : E[u,v]) : 2'E["^''l[H-*E(s)Itl 

Require: A squarefrcc polynomial / ^ 0, monic in w, and an irreducible factor e 
of discu,(/). 

Ensure: A finite set of homomorphisms E[u, i'][w] F|t] factoring through 
E[u,?;][w]/(/) s.t. the induced morphisms SpecF|i| SpecE[u,v][w]/ (f) form 
the subset of a formal dcsingularization which is centered above (e) . 
1: Fq := QF(E[u, u]/e); Let uq,vq G Fq be the natural images of u,v; 
2: if de/dv ^ then 

3: (j) : E[u, v] — > Fo[t] : u — > ug, v ^ vq + t; 
4: else 

5: (j) : Elujv] ¥o[t] : u ^ uq + t,v ^ vq; 
6: S':=0; 

7: for (cr,a) e Param(0T"'(/)) do 

8: V : Fp^i/"^)^]] -> FpJ : 7 ^ Evaluate(7, (e), t); {assuming a G F|t(i/^)^]} 
9: S := SU {E[u, v] [w] -> F[i| : u ^ ipa4){u), v ^ 'il)a(t){v), w 1-^ ilj{a)}; 
10: return S; 



In lines 1 to 5 we construct a homomorphism </> : E[m, u] Fo[t] inducing an 
isomorphism from the completed localization at (e) to Fo|t], see Corollary A. 2. 
Therefore the completions of the localizations of IC(E[u, (/)) at prime ideals 

above (e) are isomorphic to those of IC(Fo|t]/((/)^"'(/))) above (t) by Lemma A. 3. 
Lemma 3.8 tells us how these completions can be computed using a complete set 
of rational parametrizations. This is done in lines 7 to 9 using the results of a 
call to Parsun (see Algorithm 11 in Section 4.3). The homomorphism ip defined via 
Evaluate (see Algorithm 8) is just mapping t^^^ t for cosmetic reasons, i.e., 
getting rid of denominators. 

3.3.2. Divisors above Normal Crossings of the Discriminant Curve. First we fill a 
gap in Algorithm 3 and show how to test the normal crossing property for a set of 
curves. 



Algorithm 5 IsNormalCrossing(i; : 2'^["'"1) : Boolean 

Require: A set of squarefrcc polynomials E s.t. e(0, 0) = for all e G E. 
Ensure: true iff the curves defined by E are considered normal crossing at the 
origin. 

1: if E is not of the form {v, e} then 
2: return false; 
3: else 

4: return de/du{0, 0) ^ 0; 



Note that in general E describes a set of curves with normal crossing at the 
origin when E = {61,62} and det(9(ei, e2)/9(w, w))(0, 0) ^ 0. But this algorithm 
returns true only for the special situation that one of the curves is actually v = 0. 
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Remark 3.17 (Almost Minimal Jung Desingularizations). Although Algorithm 5 
does in fact test for normal crossings of a special form the overall algorithm will 
eventually terminate; Computing the point blow ups as in Algorithm 3 guarantees 
that the exceptional divisor is always of the form w = in both charts. It is known 
(see [24, Thm. 1.47]) that one can compute an embedded desingularization by a 
finite number of point blow ups. If one of the normal crossings is not of the above 
form, then our algorithm computes an additional blow up but terminates at the next 
level. In other words, the computed Jung desingularization belongs to an almost 
minimal embedded desingularization of the discriminant curve. The benefit is that 
the homomorphism in line 2 of the next algorithm is so easy to compute. One could 
also use usual normal crossings and compute more complicated homomorphisms. 



Algorithm 6 DivisorsAboveCrossing(/ : E[u,v][w],E : 2'^["'''1) : 2'E["'^lM^E(^i)[tl 

Require: A squarefree polynomial f 0, monic in w, and a set of polynomials 
E = {v,e} s.t. e(0,0) = 0, de/du{0,0) ^ and we is the squarefree part of 
disCu,(/) in the local ring E[u,u](„^„). 

Ensure: A finite set of homomorphisms E[u, F|<] factoring through 

E[u,w][u;]/(/) s.t. the induced morphisms SpecFp]] — > SpecE[M, i>][w]/(/) form 
the subset of a formal desingularization which is centered above the origin. 
1: Let g e(z,u') - u' G E[u',w'][z]; 

2: </) : E[m, v\ E|m', v'\'. V u', u t—s- ImplicitFunction((7); 
3: S 0; 

4: for (cr,a) e Param(0T«'(/)) do 

5: Let tXi, . . . , n; £ be the sequence of generators of n R>o; 

{assuming a £ E'[(u', v'Y\ and ordering generators as in Figure 2} 
6: for 1 < i < / - 2 do 

7: V : E'|(u', v'Y\ -> E'(s)M : 7 ^ Evaluate(7, (n„ n,+i), (s, f)); 

{where s,i € E'(s)pl} 
8: S := SV^ {E[u, u][w] F|t] : u ^ -^a^iu), v 1-^ 'ipa4>{v),w 1-^ 'ip{a)}; 
9: return S; 



This algorithm expects that the discriminant curve of the input surface has a nor- 
mal crossing at the origin, more precisely, the discriminant of the defining equation 
is v'^^e'^^ times a local unit. We want to compute bivariate parametrizations above 
the origin, but / is not yet normal crossing. Wc first have to apply a formal isomor- 
phism see line 2, that maps v 1— > v' and e ^ u\ because then the discriminant 
becomes v"^^u"^^ up to a unit and the defining equation 0^'"(/) is quasi-ordinary. 
To this end the image of u must fulfill the equation e{(f>{u), v') — u' = Q which has 
a unique solution by the implicit function theorem and is computed by a call to 
ImplicitFunction (see Algorithm 9 in Section 4.2). 

By the same reasoning as for Algorithm 6 the completions of the localizations of 
IC(E[m, i>][w]/(/)) at prime ideals above the origin can be computed by bivariate 
rational parametrization using a call to ParEun, see line 4. Now assume such a 
completion is given by E'[(u', d')'"]]. 

We know how to compute a formal description of a special desingularization of 
SpecE'[(M', u')'"], namely, the one obtained by a minimal number of point blow 
ups, see Lemma 3.10 and Remark 3.11. We compute the corresponding homomor- 
phisms E'[(u',w')^l E'(s)p] : {u'.v')"' ^ 5n,(m)^n,+i(m) calls to Evaluate in 
line 7. This desingularization commutes with completion by Lemma A. 4, hence, 
we get a formal desingularization of the toroidal singularity by composing these 
homomorphisms with the one given by the parametrization, compare line 8. 



18 



TOBIAS BECK 



4. Computing with Multivariate Algebraic Power Series 

In order to implement the above algorithms we have to represent multivariate, 
fractionary, algebraic power series, i.e., power series with coefBcicnts in a field E, 
variables xi, . . . ,Xn and exponents in F n R"q (the non- negative orthant of a full 
rational lattice) that are roots of non-zero polynomials in E[a;i, . . . ,x„][2:]. In the 
sequel we denote this ring by E[[x^]] (hence, this notation from now means algebraic, 
not formal power series). We must be able to perform a couple of operations: 

(1) We have to expand power series up to arbitrary order. 

(2) We have to compute the power series arising in complete sets of rational 
parametrizations for quasi-ordinary polynomials and solve equations fulfill- 
ing the conditions of the implicit function theorem. 

(3) Let hom(r,E*) denote the (commutative, multiplicative) group of homo- 
morphisms from the lattice F to the multiplicative group E* with group 
operation given by multiplication in the codomain. With cr G hom(r,E*) 
we associate the automorphism E|j:'"] ^l^Fl mapping i-^ (j{m)x"^ 
and c c for c G E. This action of hom(F,E*) on Efa;'"] should be com- 
putable. 

(4) If a G Ei|a;'"^] is a power series, ni, . . . ,n; are vectors in F/ n M"q and 
{£,1 1 ■ ■ ■ 1 a-re power series in E2[a;'"^]] where E2 is a field extension of Ei, 
we want to compute the image of a under the homomorphism : Ei — > 
Ealx^^J which maps x"^ I— > rii<i<i C"' a-iid c 1— > c for c G Ei if such is 
well-defined. 

Remark 4.1 (Origin of Requirements). Requirement (1) just means that (in order 
to return values) we need a representation in finite terms of series which a priori 
are infinite objects. 

Requirement (2) is obvious. We have to compute rational parametrizations in 
Algorithms 4 and 6 in order to compute formal prime divisors in the case of nicely 
ramified coverings. The first line of the latter algorithm also implies solving an 
equation that fulfills the conditions of the implicit function theorem. 

Further it would be nice to express certain homomorphisms Eija;'"^! E2lx'"^]. 
In general these are given by specifying the images of generators in a consis- 
tent way. But such general homomorphisms are not easy to compute. Instead we 
concentrate on the two special cases (3) and (4). 

Computing the action induced by lattice homomorphisms becomes relevant when 
computing rational parametrizations. So we need (3) to accomplish (2). 

In line 7 of Algorithm 6, we find a transformation which is a special case of 
requirement (4). Another special case of that requirement is substitution of power 
series G E2|x'"^] into a polynomial or power series g G Ei[a;i, . . . , a;„] which is an 
instance with Fi = F|j^ = Z" and the standard basis vectors. Such substitutions 
occur when we compute the composite of two homomorphisms, e.g., in line 9 of 
Algorithm 4. 

Also note that an even more special case of the last requirement is effectivity of 
addition and multiplication (with g = xi + X2 resp. g = X1X2). 

Computations with algebraic power series usually involve studying the support 
and the Newton Polygon of the defining equation (see [39, 29, 4, 35]), which live in 
Q" X Z resp. M"+^. In order to allow for a nice implementation, we will consider 
a flattened support by assuming that Q" is ordered as a group and considering 
Q" X Z a "plane" . In this setting, one can do a fair amount of "convex geometry" . 

In fact we will work with two different orderings. Let (71,(72 G Q"- We have a 
partial ordering on Q" by comparing the total degrees |gi| and \q2\ (i.e., the sum of 
their components) via < as rational numbers. Further we assume a total ordering 
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=<; which is a refinement of it, i.e., \qi\ < \q2\ imphes qi =4 q2- For a power series 
7^ Q! = J2m Q^m£™ € Els'"! we can therefore define two supports and two orders, 
namely, supp<(a) := {|m| £ Q | am ^ 0}, supp^(Q;) := {m £ Q" | am 7^ 0}, 
ord<(a) := min<(supp<(a)) and ord3j(a) := min^(supp^(a)). W.r.t. the finer 
ordering, we will also need the initial term it^(a) :— aord^{a)X°'^'^'^^"^ ■ The degree 
compatibility is important because it implies that we can easily expand series up 
to an arbitrary order w.r.t the fine ordering if and only if we can expand up to 
arbitrary total degree. 

We also have to study polynomials with power series coefficients. Consider 

0<i<d \ m / 

Its support is supp^(g) := {(m, i) e Q" x Z | g^.i ^ 0}. We define hnear maps 

(fin : Q" X Z ^ Q" : (m, i) m + in 

for all n £ Q"q and set ord^,n(g) := min^{i^n(m, z) | (m,*) G supp^(g)}. We say 
that supp^(g) has an non-trivial edge of slope n iff ipn attains its minimum for at 
least two arguments in supp^((7). Further wc define the n-th edge equation as 

edge^^n(5) := ^ 9m,rx"'z\ 

(m,?) tzsupp^ (5) and 
i^„(m,i)=ord^_„(g) 

Then supp^ (g) has a non-trivial edge of slope n iff the n-th edge equation is not a 
single term, i.e., it is also non-trivial. 

Note that supp^(g) has at most finitely many non-trivial edges. Indeed, assume 
g ~ J2i<i<k'yi^^ with 7; 7^ and 7^ ^ and let := ord^(7i) whenever 7^ ^ 0. 
Then the only possible slopes arc {i—j)~^{0j — 0i) fori > j and Oj >- O; with 7^ ^0, 
7j 7^ 0. Denoting the slopes that occur by itm we could define the "Newton Polygon" 
of g as the set {(m, i) £ Q" x I, \ I < i < k and (/3n^^ (m) ^ ord^_ri,„ (<?) for all m}. 

4.1. Representing Algebraic Power Series. In the sequel we suggest a repre- 
sentation which is suitable for the computer algebra system MAGMA [7] and facili- 
tates all of the above operations. We only highlight the essential points. MAGMA, 
like most other computer algebra systems, does not provide polynomials / e E[a;^] 
with fractionary exponents directly. For simplicity of reading, we nevertheless use 
such polynomials, an implementation being straightforward. We also assume that 
given / we can ask for its coefficient field E and exponent lattice F and that an 
implicit conversion mechanism is provided for Ei[x'"i] C E2[a;^^] whenever Ei C E2 
and Fi C Further it should be understood that a tuple defining a series may 
recursively depend on other series or even on polynomials with series coefficients. 
We represent algebraic power series using a hybrid lazy-exact approach by finite, 
acyclic, directed and rooted graphs with nodes of two types: 

Type A: We represent an algebraic power series a G E2|x'"^| by a tuple 
(4.1) (ao,/) 

where ao G E2[x^^] is an initial segment of a w.r.t. =^ and / = J2o<i<d 'fii^^ ^ 
lEi[x'"^][2;] are such that Ei C E2, Fi C r2 and / ^ is squarefree (when 
considered as a polynomial over QF(Ei|£^i])) and vanishing on a. 

Such a node in general has as c?-|- 1 descendants, namely, the coefficients 
(fii. As a special case we allow / £ Ei[a;^i][z]. In this case we store / itself, 
there are no descendants and the node is terminal. 
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We also need a technical condition to ensure that oq identifies a uniquely 
amongst the roots of /. 

Algorithm 7 Series(ao : E2[x^'']J ■ Eife^iW) : E2k^'l 

Require: A tuple as in (4.1) fulfilling Condition 4.3. 
Ensure: The power series a defined by it. 
1: Encapsulate data and return object of type "algebraic power series"; 



Type B: The second type of node is given by a tuple 

(4.2) (a,n,C) 

and represents the image (3 :— 4>{a) under the homomorphism described in 
the third requirement. It has 1 + 1 descendants, namely, the and a. Also 
in this case we need a technical condition necessary for (f> being well-defined. 

Algorithm 8 Evaluate(a : Ei|x^i],n : (F/ nR!^o)',C : Ezfe^^J) : Ealx^^J 

Require: A tuple as in (4.2) fulfilling Condition 4.6. 
Ensure: The power series f3 defined by it. 
1: Encapsulate data and return object of type "algebraic power series" ; 



Nodes of type B facilitate explicitely the operation of requirement (4) from our 
list. The following sections describe how to implement the remaining requirements 
using algorithms based on an idea in [3] which is similar in spirit to [35] . 

Remark 4.2 (Representation Paradigm). The reason for calling this representa- 
tion lazy- exact is the following: In order to expand a series up to some order, we 
will recursively expand its descendants (maybe storing precomputed values) and 
then compute the approximation of the series under consideration in a lazy fash- 
ion. On the other hand, using elimination theory and traversing the representation 
graph recursively, one can compute a (minimal) defining polynomial for any series 
represented as above. Therefore zero-equality and polynomiality /rationality are 
decidable. In this sense we speak of an exact representation. Since in particular all 
polynomials (including and 1) are representable we get a computationally effective 
ring. All this (and a little bit more) is included in the MAGMA implementation. Note, 
however, that it is advisable to avoid the effective decision algorithms because they 
depend on nested resultant computations and can be very expensive. 

4.2. Expanding Algebraic Power Series. Now we show how to expand power 
series in a lazy fashion. For more efficient computations with lazy power series, see 
[34]. For expanding a power series a given by a tuple (4.1), we consider 

d 

(4.3) g := f{z + uq) ~ ai := a — ao, n := ord^(Q!i) and c := ord3;((7i). 

i=0 

First we want to find implications of the fact that g{ai) = 0. The case ai = 
happens if and only if 17(0) = or, equivalently, 70 = 0. Let's assume ai ^ and 
set a' := it=<;(Q;i). Next we study the contribution of the terms 7^2* to the result 
under the substitution z 1— > ai. Whenever 7^ 7^ we find 

jia\ = 7i(a' + (cvi — a'))' = it^(7i)a'* + . . . (higher order terms). 

The minimal order of these expressions is 

min{ord^(it^(7i)a'*) \ < i < d and ji 7^ 0} = ord^^n(5)- 

ord^ (7i)+in 
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Since g{ai) = and a' 7^ it follows that the terms of order ord^.n(5) must sum 
up to zero. In particular, there have to be at least two indices i s.t. the terms 
it^(7i)a'* are of order ordac;,n(<?)- In other words, the n-th edge equation of g must 
be non-trivial and a' must be a root of it. 

If we can make sure that the data in (4.1) implies that the n-th edge is linear 
of the form it-j(7i)z + it^(7o) = then a' is uniquely determined and easy to 
compute. 

Condition 4.3 (Valid Representations of Type A). With definitions as in (4.3) 
we require for valid representations of type A that either 70 = or edge^ „(g) is 
linear. 

Now substitute ao ^ ao + ct' , ai ^ ai — a' and g ^ g{z + a'). A careful 
analysis of the exponent structure shows that under the above condition it^(7i) 
doesn't change and Condition 4.3 remains valid. Hence we can repeat the argument 
and compute successively as many terms as needed. 

Example 4.4. Let f :— — Zxiz^ — ^x\x\z^ -V "ix^z^ — x\. Then / has a power 

1/2 1 2 /3 

scries root a starting with ao = x^ + "^x^ X2- We work with the degree graded 
lexicographical ordering =:<;. First we translate / as in (4.3): 

f{z + ao) =z^ + + lxl''^X2)z^+ 

{+12X2 - 'ix^xT + ixrxj)z^+ 

( — 8X2^ + %X^ X2 ~ 15^1^ xj. + '^2&'^^^'i^'^ ~^ 

OXj^ X2 -|- gX]^ X2 4096 1 X2jZ -^ 

4/3 7/2 J_ 2 4 , _9_ 8/3 9/2 9_^10/3^5\^ , 

V 8 1 2 "T 54X^X2 -r 2048 1 2 16384 1 ■^2)'^'^ 

I 3_ 8/3 5 , _9_ 10/3 11/2 ^^4 6n 

V 1024-^1 -^2 "T 16384"^1 "^2 262144 '^l'^2j 

We find that /(z-|-ao) has a hnear edge of slope (|, |) and extract the edge equation 
— l^^^'^Xj^^z — Tc^^i^"^^! — 0- The solution is the next term — j^gaj^^'^Xj^^ in the 
expansion of a. Next we consider f{z -I- ao — Y2»x\'^ :^2'^') and find that it has a 
linear edge of slope (|, |). Solving the edge equation —\x^^'^X2'^ + 262^44 ^1^2 = 
we find the next term 32768 ^1'^ '^-'^2^^' ^"^^ on: 

_ 1/2 1 2/3 \_ 4/3 3/2 1 8/3 5/2 _ 1 4 7/2 

a — X2 -Y ^X^ X2 128 '''1 ^1 ' 32768 '''1 '^2 4194304 •''1"^2 + ■ • ■ 

Obviously we do not need to know g completely to do this computation. More 
precisely, we can expand a up to order less than o, if we have approximated g 
sufficiently well. All the terms of order less than o arc determined by a linear edge 
equation of the form \X^{^\)z + ■ ■ ■ = 0. Therefore, it is sufficient to know the 
constant term up to order less than ord3;;.(,(it^(7i)z) = c + o with c as in (4.3). 
Hence, for approximating g it is sufficient to expand the coefficients of / up to 
order less than c -f o. 

For example. Condition 4.3 is true when considering a power series defined by 
the implicit function theorem. Therefore solving such equations is trivial: 



Algorithm 9 ImplicitFunction(5 


: Eilx^^i 


lW):]Eik^1 


Require: A polynomial g s.t. 17(0, . . 


.,0,0) = 


= Oand ag/az(o,...,o,o)^o. 


Ensure: The unique root a of g s.t. 


a(0,... 


,0)=0. 


1: return Series(0, 5); 
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Indeed, if g is equal to J2i=o li^^i 3(0; ■ • ■ , 0, 0) = and dg/dz{0, . . . , 0, 0) ^ 
then ord^(7o) >- (or 70 = 0) and ord^(7i) = 0. If 70 ^ then edge^ ,,(5) must 
be Uncar where n — ord3^(7o) = ord3^(a). 

Remark 4.5 (Universality of Type A). Note that for any algebraic power series 
a there is, by definition, a squarefree polynomial / with polynomial coefficients 
vanishing on it. This implies that f{z + a) vanishes at z = with multiplicity one 
or that the constant coefficient of f{z + a) is zero whereas its linear coefficient does 
not vanish. This means that the initial term of the linear coefficient is fixed if we 
consider translations f{z + ao) by sufficiently large initial segments ao of a. As a 
consequence we can always find an initial segment ao s.t. Condition 4.3 is fulfilled 
and every algebraic power series is representable by a single node of type A. 

So our representation is already universal. But computing efficiently the last op- 
eration of our requirement list is not easy in this representation. So for algorithmic 
purposes we also allow the second type of node. 

The previous discussion showed how to expand power series represented by a 
graph consisting of nodes of type A. Since =<; is compatible with total degree this 
implies that we can expand such series up to arbitrary total degree. Now we 
consider a power series /3 represented by a tuple (4.2). Indeed what comes next is 
best explained in the total degree ordering. 

Assume we want to expand (3 up to order less than o G Q. Therefore we first 
compute approximations of the descendants up to order less than o. Then 
we apply the map ni<i<; ^"'^""^ to each term in an expansion of a, sum up 

the intermediate results and finally truncate at order o. For this truncation to be 
correct, we have to use a sufficiently large expansion of a, say, up to order less than 
o' . The only remaining question is how to determine o' from o. 

To this end let 

(4.4) " E °^'^<(^«) = ("1' ■ ■ ■ a-nd c ;= min{nj | 1 < j < n} 

\i<'<' / 

where ord<(0) := 00. Let </) denote the homomorphism x"^ i~> ni<i</ C^"^^ then 



This calculation shows two things: First, if c > then for any 7 ~ X^m^'TiS"' the 
sum 7m0(i£") converges, so (t>{^) is well defined. This is an analogue to the 
usual condition when substituting formal power series into each other. 

Condition 4.6 (Valid Representations of Type B). For valid representations of 
type B we require c > for c defined as in (4.4). 

Second, under this condition, if a = J2m '^mSl^ then the terms amiS™ with |m| > 
o/c contribute terms to (j){a) of order greater or equal o, so we can choose o' := o/c. 

Remark 4.7 (Contraction Constants). We can give a nice theoretical meaning to the 
two technical conditions and the deduced algorithms. Namely, the values c and c 
may be understood as additive resp. multiplicative contraction constants for certain 
continuous maps between power series domains with the usual metrics. In the 
case of Condition 4.1 this map is contractive only in a small enough neighborhood 




n(m) 




cmj = c\m\. 



l<j<n 
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(determined by ord^(a — ao)) of the root a. Contractivity makes the represented 
series well-defined. 

In both cases the constants determine how far the descendants of a node have to 
be expanded. An implementation could store them together with the representing 
tuples (4.1) and (4.2). We would also like to mention that for an efficient imple- 
mentation it is crucial to truncate intermediate results (e.g. defining polynomials or 
approximations of descendants) as often as possible to a sufficient precision. The 
bounds are computed by similar calculations as above. We omit them from the 
presentation and show only the essentials of the algorithm. 

4.3. Rational Parametrizations for Quasi-Ordinary Polynomials. A quite 
intricate thing is the computation of complete sets of rational parametrizations. 
In requirement (3) we have introduced the action of the group hom(r,E*) on 
E|x^]]. We write this as a left action using again functional notation, i.e., for 
a G hom(r,E*) and a G Els'"] wc simply write cr{a). This is an exponent structure 
preserving automorphism in the sense that supp^(Q:) = supp^ (cr(Q!)). Let's con- 
vince ourselves that this action can effectively be carried out in our representation. 

Remark 4.8 (Effectivity of Actions Induced by Lattice Homomorphisms) . Since 
such an automorphism is structure preserving, it is easily computed in our rep- 
resentation: If a power series is represented by a tuple (4.1), then we apply the 
automorphism to the initial segment ao and to the coefficients of the defining equa- 
tion /, possibly recursing to the descendants. If a power series is represented by a 
tuple (4.2), wc apply the automorphism recursively to the ^i. 

Definition 4.9 (Rational Parametrizations). We call (cr, a) with a G Ei|x^i] 
and a € hom(ro,E*) a parametrization of a monic polynomial / G Eq [a;'"" | [z] iff 
Eq C El, Fq C Fi and a^^'{f){a) = 0. The order of the parametrization is defined 
to be ord^(a). 

Let g\f be an irreducible factor s.t. tT^^(f/)(a) = 0. We call (cr, a) rational if 
the induced homomorphism IC(Eo|x'"''][2;]/ {g)) 'Eifa^'^l which maps z i— > a and 
7 (7(7) for 7 G EqIx'""]] is an isomorphism. 

Intuitively, rational parametrizations are distinguished by minimal field and lat- 
tice extensions. The induced homomorphism exists due to the universal property 
of integrally closed domains. If this homomorphism is an isomorphism then the 
irreducible polynomial g from above at least has to be prime. (Mind that EqIx'""] 
is in general no UFD!) 

Remark 4.10 (Irreducible Monic Polynomials). Let g G Eo|x'"°|[z] be a monic 
irreducible polynomial. We want to prove that g is prime. The following arguments 
are taken from [28]. It is sufficient to show primality in QF(Eo[a;'"'']])[2;]; Indeed, 
if hg G Eo|x^°][z] for some h G QF(Eo[i'"''|)[2:] then a Gaussian style inductive 
argument shows that h G Eo|x^i[z]. So g EqIx^^Hz] = (5 QF(Eo|xi^»l)[z]) n 
Eo|x'"°][z], i.e., it is the preimage of a prime ideal and therefore prime itself. 

Let g' G QF(Eo[x'"''])[z] be an irreducible (and hence prime) monic polynomial 
factor of g. We show g' = g. Being also roots of g, all roots of g' (in some 
splitting field) are integral over EqI^'"'^]. The coefficients of g' (being polynomials 
in these roots) are also integral over EqIx'""] and elements of QF(Eo[x'"'']). Then 
g' G Eo|x'""][z] because Eo|a;'"o| is integrally closed. By the above argument g'\g 
also in Eg [x'"'' ] [z] and so g' = g because g is irreducible. 

This also implies that a monic polynomial / G Eo|iE'"°][2:] has a unique (up to 
permutation) factorization into monic irreducibles. 

Let now / G Eo|x'"°][z] be quasi-ordinary and assume we want to compute a 
complete set of rational parametrizations. Let's be more general and say we want 



24 



TOBIAS BECK 



to compute a complete set of rational parametrizations of order greater than some 
value no G Q". (With the choice no ( — 1, . . . , —1) this includes the original task.) 

First assume we are given a parametrization (ct, a) of / s.t. a ^ and n := 
ord=^(Q;) >- no. Write g := cr^^{f) and a' := it=<;(Q;). By the same reasoning as in 
the beginning of Section 4.2 (with a in place of ai), wc find that the n-th edge 
equation of g must be non-trivial and a' must be a root of it. 

Now assume we have to find (cr, a) as above using a recursive approach. As a first 
step we have to determine the initial term of a (up to some isomorphism). From 
the previous discussion it follows that its slope will be n ;^ no s.t. / has a non-trivial 
n-th edge equation. If the tuple (cr, a) is meant to be a rational parametrization, 
then we have to solve this equation economically, i.e., using a field extension of 
least degree, and at the same time determine a partially. Duval's trick [12] adapted 
to the multivariate case gives an optimal choice. 

To this end let for 1 < i < n be a basis of Fo. Let b E Z>o be the size of 
F'/Fo where F' := Fq + Zn, i.e., b is minimal s.t. 5n G Fo. Now we can write 

edge^,,(/)= ^ /,x"'z' = x"'z' fi+Jb^-''"^" 

l<i<k 0<j<(k~l)/b 

with fi ^ and fk ^ 0. Further we can express bn — X/i<i<n '^i'^'i where nec- 
essarily gcd(fo, ci, . . . , c„) = 1 since b was chosen minimally, and compute Bczout 
coefficients u and Vi s.t. u^+^K^n ^i^i — 1- (Note that this choice is not unique.) 
If r is a non-zero root in a minimal field extension E' s.t. 

E /'+^'''^' = 

0<i<(*:-0/'' 

then we can define the homomorphism cr' : Fq E'* : ^ r~'"' and the initial 
term a' := r"x". With these definitions one checks 

edge^,„(a'T^(/))(a') = a'T^(edgc^,„(/))(r"x") 

0<j<{k-l)/b 
^<j<{k-l)/b 
n<j<{k-l)/b 

= (•••) E f^+l" v.c^-bn+bny 

0<j<{k-l)/b 

^(...) E W-'=o. 

0<j<{k-l)/b 

This finishes the description of how to choose the first term (up to isomorphism) 
and part of the structure preserving automorphism. The remainder of course has to 
be chosen via recursion. More precisely, define /' := a''^^ {f){z + a') and compute a 
parametrization (cr", a") S hom(F',E*)xEi|x^i]] of /' with ord^(a") >- n. Then we 
may just set a := cr"(ci;') 4- a" G Ei|x^i| and a := a" a' G hom(Fo,E^). (Note that 
we have an inclusion hom(Fo,E'*) — > hom(Fo,E*) and a surjection hom(F',E*) 
hom(Fo,EJ), hence, we can build a"a' G hom(Fo,E^).) Now trivially 

= a"T^(/')(a") = a"T-(a'T-(/)(z + a')){a") 
= f7"T^(a'T-(/))(z + <j"{a'))ia") = a"T-(a'T-(/))(a"(a') + a") = a^^{f){a). 
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Example 4.11. Let /o := + 3x2z'^ + x\x2Z^ + 2>x^z^ + xf. Its discriminant is 
disc4/o) = 729xf 64) and, hence, it is quasi-ordinary. We again use the 

degree graded lexicographical ordering =<;. 

The only rational slope of a non-trivial edge is ni := (0, ^) with edge equation 
-\-'ix2z'^ ^ix^z^ + x^ = 0. We sec that only even powers of z have a non- vanishing 
coefficient which corresponds to the fact that 2ni e 1? , i.e., hi = 2. Let ri be a 
solution of = + 'ir^ + 'ir + l = {r + lY , hence, ri = —1. The standard lattice has 
basis ~ (1, 0) and mi, 2 = (0, 1). Then 2ni,i = Omi,i + lmi,2, hence, ci,i = and 
C1.2 = 1- One verifies the Bezout relation O61 -|-Oci.i + lci,2 = 1 and we have ui = 0, 
= and wi,2 = 1- Therefore we define the first term ai := r^^Xj^^ = 12^^ 
and the homomorphism to be CTi : mi,i 1— > "^'^ = l,Tni_2 ^ '^1 "^'^ = —1 with 
corresponding action xi 1-^ Xi,X2 '-^ —X2- We set 

/i:=^i^^(/o)(;^ + 4) = 

1 3 7 8 9 

z^ + 6x1 + I2x2z^ + {8x^ - xlxl)z^ - Sxjx^ z^ - ?,x\x^ z - x\x^ . 

For /i we now find two rational slopes of non-trivial edges, namely, rii from above 
and n2 = (§7 !)• The extended exponent lattice at this point is Zm2,i +ZTn2.2 where 
tn2,i = (1, 0) and m2,2 = (0, ^). We have 3n2 = 2m2,i -I- 6m2,2, i-e., 62 = 3, C2,i = 2 
and C2,2 = 6. From a Bezout relation we get U2 = 1, 1^2,1 = —1 and W2,2 = 0. The 
edge equation is 8x2^^ z'^ ~xfx2^^ = 0. Let r2 be the root of 8r — 1 = 0, i.e., 7-2 = |. 
Hence, we define a new homomorphism cr2 : m2,i 1-^ "^'^ = ^,^2,2 f2^'^'^ = 1 
with action xi 1— > ^xi,X2^^ 1— > X2^^. We also define a new initial segment a2 '■= 

o-2(q;i) -I- r2^xJ/^X2 = x^^^ - ^xl^^X2 and set /2 := C72^''ifi){z - ^xl^^X2). 

From now on the edge equation will always be linear and we have determined 
a parametrization. More precisely, the overall automorphism is cr2f 1 and acts by 
xi 1-^ \xi,X2 I— > —X2. The power series starts with 02 and is a root of cr2cri^^(/o)- 
We have seen its expansion in Example 4.4. 

This discussion yields the following algorithm. The return type Eo|x''^ ] means 
algebraic power series with coefficients in some finite algebraic extension Eg and 
exponents in some rational lattice containing Z". 



Algorithm 10 ParamRec(/ : Eolx^^i [z], no G Q") : 2'^°™(ro,Eo*)xEoy;'0"l 

Require: A quasi-ordinary polynomial / and an order Hq. 
Ensure: A set of rational parametrizations of / of order greater no. 

1: 5* := {n G Q>o | n >- no and edge^ ^{f) is non-trivial}; V ■= 0; 

2; if (Vn e 5 : ord,(edge^ „(/)) > 1) or (3n G 5 : dcg,(cdge^ „(/)) = 1) then 

3; 5 5 \ {n G 5 I deg,(edge^ _„(/)) = 1}-V:=V\J {(1, Series(0, /))}; 

4: for n G 5 do 

5: Let b G Z>o be minimal s.t. bn G Fo; Let for 1 < j < n be a basis of Fo; 
6: Write bn = X]i<i<n '^i^i ^-^i^ compute u,Vi G Z s.t. ub + J2i<i<n '^i'^i — 
7; Write edge^ „(/) = x^' z' Eo<,<fc f,3l'""'z^'' with fo + and ^ 0; 
8: R := ZeroSet((^o<^.<^ f^z^)); 
9: for r G i? do 

10: a' := r^x"^; cr' : Fo ^ E'* : ^ r^^'; {assuming r G E'*} 

11: P:=VU map(^„_„„)„(^„^, ,,„(„,)+„„)ParamRec(cr'f^(/)(2 + a'),n); 

12: return V: 
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For the definition of the auxiliary function ZeroSet see Remark 3.15. Applying 
structure preserving automorphisms and translating the defining polynomial in z 
preserves quasi-ordinariness. Therefore the arguments to the recursive call always 
fit to the input description. Although the algorithm never makes explicit use of 
Theorem A. 5 (the Theorem of Jung-Abhyankar) it depends crucially on / being 
quasi-ordinary. For example, the power series constructed in line 3 are well-defined 
in that case; Otherwise there might not be any power series root of / which is 
supported on the positive orthant and whose initial term is a root of the linear 
edge of /. Recursion in this algorithm ends when the set S is empty when entering 
line 4. Termination is assured by the very argument of Remark 4.5. 

Remark 4.12 (And a Little Bit of Engineering. . . ). In line 1 of the above algorithm 
we compute the "Newton Polygon" of / = J^i fi^^- Since / is monic, one of the 
non-trivial edges contains the point ((0, . . . , 0), deg2(/)). Now z divides / either 
with multiplicity or 1 (because of squarefreeness). So for computing the non- 
trivial edges it will be enough to expand the coefficients ipi until the initial term 
of ifi or ipo appears depending on whether z divides / or not. In other words we 
would have to check, whether ipo = 0. In principle this is possible, sec Remark 4.2, 
but from the point of efficiency it is not advisable. 

On the other hand it doesn't matter if z divides / or if / has a linear edge 
equation (compare to the condition in line 2). In both cases we return Series(0, /) 
which might be zero. So an engineering solution might be to compute either all 
non-trivial edges or only a set of non-trivial edges with lowest vertex at linear level 
and make sure that the constant term (if existent) has large enough order. 

Another approach would be to modify the above algorithm to allow approximate 
input with polynomial coefficients only but also include error reporting when an 
approximation was not accurate enough. 

Further it is advisable to choose a short vector (u, wi, . . . , w„) G Z"+^ using for 
example the LLL-algorithm. For practical purposes the additional complexity is 
negligible compared to the grow of coefficients if a larger vector is chosen. 

We didn't include the technical tricks into the description because the algorithm 
is easier to read and to argue about when written up as above. 

Setting now no := (—1, . • . , —1) we get a set of parametrizations of /: 

Algorithm 11 Param(/ : Eo[£ro][z]) : 2h°m(ro,EF*)xE^Iir«"l 

Require: A quasi-ordinary polynomial /. 
Ensure: A complete set of rational parametrizations of /. 
1: return ParEiinRec(/, (—1, . . . , —1)); 



We have to show that the computed parametrizations are rational and complete. 

Lemma 4.13 (Bounding Extensions). Let Algorithm 10 be called with a quasi- 
ordinary polynomial f G Eq ]] [2] and no G Q" and write {ui, ai) S hom(ro, E*) x 
Ei|x^'] for the computed parametrizations. Then we may hound the extensions from 
above by 

^[E, : Eo] #(r,/ro) < max{deg,(edge^ _„(/)) | n >■ no}. 

i 

Proof. Assume the condition in line 2 holds. Then cither / has a zero root, i.e., 
z\f or we have found a linear edge equation. In both cases we construct one 
parametrization of / as in line 3. The involved series is an element of Eo[a;'"''] and 
the lattice remains unchanged. This parametrization contributes a summand 1 to 
the lefthand side of the inequality. 
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Now set d := niax{deg_j(edge^ n(/)) I ^ ^ "^0} — 1 if the condition in line 2 was 
true, d := max{deg^(edgc^ n(/)) I ^ ^ ^0} otherwise. Wc have to ensure that the 
sum modified by running only over the parametrizations constructed in the loop is 
bounded from above by d. To this end assume that S contains o slopes nj when 
entering line 4. Let £j denote the height of the edge equation of slope xij, i.e., the 
difference of its z-degree and its 2-order. Then J2i<j<o — ^■ 

Fix an edge, i.e., a value of j. Then we compute an integer bj in line 5 and a 
number of roots rj.fe S ^j,k with multiplicities, say, ruj^k for 1 < k < Sj in line 8. 
It follows that bj — ^(Tj/To) with Tj := Fq + Zn^ measures the extension of the 
exponent lattice and £j = bj{J2i<k<s l^J-k ■ ^0] 

Fix a root, i.e., a value of k and let fj^k be the polynomial used as parameter 
in the recursive call in line 11. Now rrij^k is also the multiplicity of the root of 
the corresponding edge equation. After translation zero becomes an m^-^fc-fold root, 
hence, the 2-order of edge^ is equal to mj_k. This can serve as an upper 

bound in the statement for the recursive call, which returns, say, tj,k different 
parametrizations . 

Assuming that the lemma is true for recursive calls (sec (*) below) one computes: 

l<j"<o l<J<o l<fe<Sj 

l<j<o l<k<Sj i<l<tj.k 

= J2 [E,-,,i :E,,,][E,,fe:Eo] #(F,,fe,/F,)#(F,/Fo) 

i<j<o i</£<,sj i<;<tj fc 

= E E E [E,,M : Eo] #(F,,fe,/Fo) 

i<j"<o i</£<,sj i<;<tj,fc 

For each tuple of indices {j, k, I) appearing in that sum, the loop now produces 
exactly one parametrization with coefficient field Ej^k,i and exponent lattice Tj^k,i- 

' '□ 

For what follows we assume that Eq C Ej C Eq. A parametrization {ai,ai) € 
hom(ro,E*) X Ei|x^'| is very close to a root of / in Eo|a;^'|. It provided a root 
if we could reverse the effects of cr^. Precisely, let r G hom(ri,Eo) be an extension 
of the inverse of Ci to F^, in other words, tItq = <^i~^ ■ Then we call /? := T{ai) 
an associated root to (17^,0!^), because applying r to the equation = aj^(/){ai) 
yields = {Ta^y''{f){T{ai)) = ((Ti'V,)T^(/)(/3) = /(/3). It is not hard to show that 
can always be extended to hom(Fi,Eo) because Eq is algebraically closed. 

The next lemma says, that we do not miss any associated roots. 

Lemma 4.14 (Completeness). Let Algorithm 10 be called with a quasi- ordinary 
polynomial f G Eo|x'"''|[z] and rio G Q" and write {(Ti,ai) G hom(Fo,E*) x 
Eilx'"'] for the computed parametrizations. Then each root (3 G Eo|x^ J of f 
with ord^(/3) >- Wq is associated to at least one (ai,ai). 

Proof. Let f3 he a root of / of order n :^ Hq. If /3 = or if /3 ^ and edge^ „(/) is 
a linear equation then we consider the parametrization constructed in line 3: The 
involved series is already f3 and the involved automorphism is the identity. 

Otherwise assume the outer loop is processing slope n, set F' := Fq + Zn and 
write /? = /?' + (3" where /3' := it=<;(/3). Then f3' must be a root of edge^_„(/). 
Setting /?' = ox" for some o G Eq we have that r :— o^ is a root of X]o<i<fc /j-^"'' 
compare line 8. 
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Assume now we are in the inner loop processing this r G E'* and let /' := 
a'^^{f){z + a') be the polynomial in the arguments to the recursive call in line 11. 
Let t' € hom(r',Eo) be defined via t'Itq = cr'"^ and n r~'^o. We have to check 
that this is well-defined, more precisely, since bx\ S Fq (with b minimal) we must 
have that (r~"o)'' coincides with (T'^^(6n); Indeed 

In particular T'{a') = r'(r"a;") = ox" = /3' and applying r'^^ to = /(/3) we get 

= (r'-i)T-(/)(r'-i(/3') + r'^n/?")) = +^'"1 = /'(r'-^ (/3"))- So 

r' ^(/3") is a root of /' of order greater n. 

Assuming that the statement holds for the recursive call, we get a parametriza- 
tion (cr", a") € hom(ro,E"*) x E"|xr"]l which is associated to T'-\j3") via, say, r" 
with r"|r' ~ a"~^ . This is combined to a returned parametrization (ct'V, ct"(q;') + 
a"). We claim that this parametrization is associated to /? via t't" . Indeed t't" 
restricts to the inverse of a" a': 

(r'r")|r„ = (r'|ro)(^"|r„) = o'-\o"-\,:) = {a" a'Y' 

It remains to show that it maps a" {a') + o" to /3: 

r'T"(a"(a') + a") ^ T'(r"a")(a') + r'(r"(a")) = 

r' (a"- 1 a") (a' ) + r' (r'" ^ (/?" ) ) = r' (a') + /?" = /?' + /3" = /3 □ 

Now if /3i and /32 are roots of / associated to (cr,;, Ui) via ti and T2, then t\t^^ 
restricts to the identity in hom(ro,Eo), hence, maps Fq to 1. In other words t\t^^ 
acts as an automorphism on Eo|x'"']l fixing Eo|a;'"'']. Therefore, the roots which are 
associated to one and the same parametrization are all conjugate and, hence, roots 
of the same irreducible factor. This gives an injection from the irreducible factors 
(with roots of order greater no) to parametrizations (of order greater no). 

Corollary 4.15 (Complete Sets of Rational Parametrizations). Let Algorithm 11 
be called with a quasi- ordinary polynomial f G Eq [i'"" ] [z] and write {ai,ai) G 
hom(Fo,E*) X Ei|x^'] for the computed parametrizations, then: 

• deg,(/) = EJE, :Eo] #(F,/ro) 

• All computed parametrizations are rational and in bijective correspondence 
with the irreducible factors of f . 

Proof. To show this, we apply the above two lemmata with no := (— 1, . . . , — 1) 
each. In Lemma 4.13, clearly, dcg^(f) is an upper bound for the righthand side of 
the inequahty, so degj/) > EJE, : Eo] #(r,/Fo). 

Now let fj\f be the irreducible factors. By Lemma 4.14 and the previous dis- 
cussion we may assume that there is an injection j i-^ ij s.t. fj is parametrized by 
(cTi^- , ) . Then we have chains of field inclusions 

QF(EoIair°]) C qF{Eolx^"M/ if,)) C QF(E,j/'.l), 

where the second inclusion is given by (ai . ,ai). The extensions are algebraic of de- 
grees [QF{Eolx^mz]/ {f,)) : QF(Eo[a;r"])] = degj/,) respectively [QF(E,Jx^'. 1) : 
QF(Eo[£r"])] = [E,^. : Eo] #(F,7Fo) and therefore 

deg,(/) = ^deg,(/,) < ^[E,^ : Eo] #(F,^/Fo) < ^[E, : Eo] #(F,/Fo). 

j j i 

Together we have proven equality. In fact, since all summands are positive we 
find that the map j i-^ ij is a bijection and for each chain of inclusions as above 
equality of degrees must hold. This gives the second statement. □ 
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5. Conclusion 

We have introduced the concept of formal desingularizations and shown how 
to compute them for hypersurfaces of Pj|. The algorithm has been implemented 
and found to run very well. We can compute formal desingularizations faster by 
magnitudes than the general algorithms for dcsingularization. The reason is that 
our algorithm doesn't depend on Groebner basis computations. It relies only on 
linear algebra and polynomial factorization. 

Here, we would like to point out an application of our algorithm. Namely, formal 
desingularizations can be used to compute adjoint sheaves. We will describe how to 
do that in another paper. Adjoint sheaves in arbitrary dimension arc defined to be 
the direct image of the tensor powers of the canonic sheaf w.r.t. an arbitrary dcsin- 
gularization. Adjoint spaces for projective schemes (i.e., graded components of the 
associated graded ring to that sheaf) are the keystone to the rational paramctriza- 
tion of curves. In [32] it has been shown that they are of equal importance for 
the computation of rational surface parametrizations. For example, they facilitate 
the computation of the arithmetic genus and the plurigenera of the surface. Thus, 
we can effectively check Castelnuovo's Criterion for the parametrizability of sur- 
faces. Moreover, adjoint spaces can be used to construct certain rational maps 
that reduce the paramctrization problem to a set of base cases. The final goal is 
an efficient implementation to rationally parametrize hypersurface in Pq (with or 
without introducing field extensions). 

We finish with an open problem. With our current definition, a formal dcsingu- 
larization is just a loosely related set of formal prime divisors. For some applica- 
tions, however, it would be nice to know the dual graph of the surface dcsingular- 
ization, i.e., an annotated graph with one vertex for each exceptional divisor and 
edges whenever two divisors intersect. Such a graph could be the starting point for 
an algorithm to compute the minimal formal dcsingularization. In our method it 
would be easy to keep track of the dual graph of the embedded curve dcsingular- 
ization and in a certain sense the dual graph of the resolution projects to it. Such 
graph coverings have been studied in [30, 31]. The problem is that we get formal 
prime divisors from two sources, namely, Algorithm 4 and Algorithm 6. However, 
it is not clear how the corresponding prime divisors intersect. 

The algorithm of this paper as well as the method for the computation of adjoint 
spaces is available as a Magma-package and can be downloaded via this link: 

http : //www. ri cam. oeaw. ac . at/sof tware/symcomp/adjoints . tar .gz 
It will probably become part of the next major Magma-release. 
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Appendix A. Some Local Algebra 

In this appendix we gather a few results from local commutative algebra and 
present them in a form suitable for our needs. Again E denotes a field of character- 
istic zero which need not be algebraically closed. We are dealing with completions 
of stalks of regular schemes of finite type over E. Therefore we first give a famous 
structure theorem in this setting. Recall that an essentially finite local E-algebra 
is the localization of a finitely generated E-algebra at a prime ideal. 

Theorem A.l (Cohen Structure Theorem). Let {A,m) be the completion of an 
essentially finite, regular, local algebra of Krull dimension s. SetW^ := A/m with 
canonic projection tt : ^ — > Fq. Further let 
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• {ui}i<i<r a A be a set projecting to a transcendence basis o/Wq over E and 

• {^j}i<i<.5 C m a minimal set of generators. 

Then there is a unique coefficient field F C yl containing E and {ui}i<i<r s.t. n 
restricts to an isomorphism F ^ Fq and 



F|xi,...,a:J 




for few, 

for I < j < s 



is also an isomorphism. 



Proof. The existence of a unique field F fulfilling the first assertion is the content 
of [27, Thm. 28.3] and its proof (where references to differential bases can be sub- 
stituted by transcendence bases in characteristic zero). Since {wj}i<j<s is a set of 
generators for m one easily sees that the homomorphisni in the second assertion is 
surjective, cf. [2, Lem. 10.23]. The rings on both sides have the same dimension. 
Therefore the kernel must be trivial by [2, Cor. 11.18]. □ 

We actually need a constructive version of a kind of inverse of the above isomor- 
phism in a special case. 

Corollary A. 2 (Completion at Points in AfRne n-Space). Let p := (/i, . . . , fr) C 
E[.Ti, . . . , x„] be a prime ideal of height I < min(r, n), 



rr 9{fl, ■ • ■ , /r) ^ 



r X n 



d{xi, 

the Jacobian matrix, ¥q := QF(E[a;i, . . . ,Xn]/p) the residue field, tt : E[a;i, . . . ,Xn] — > 
Fq the canonic projection and 

(/) : E[xi, . . .,Xn] Fo[ti, ...,ti] : ^ Tr{xi) + Mi (<i, . . .,ti)'^, 

where Mi are row vectors with entries in {0, 1}. Write M := (Mi)i<i<„. 

Then 4> extends uniquely to a homomorphism E[a;i, . . . , a::„]p — > Fo|ti, . . . , t;]. 
Moreover we can choose M s.t. ■n{J') M £ Fq'^' has rank I and then the extended 
homomorphism becomes an isomorphism. 

Proof. The homomorphism trivially extends to Fo|ti, . . . ,t;]]. Then for any g G 
E[a;i, . . . , Xn] one computes for the image 

0(5) € T:{g) + T:{dgld{xi,...,Xn))M{ti,...,tif + 

where m := {ti, . . . ,ti). Now the above equation shows that if g e E[a;i, . . . , x„] \ p 
then we have TT{g) ^ for the constant part. Therefore the image is a unit and 
the homomorphism lifts uniquely to E[a;i, . . . ,x„]p. Impose the p-adic topology 
on the domain and the m-adic topology on the codomain. Writing down the long 
expansion one sees that the homomorphism is even uniformly continuous. Since 

the codomain is already complete we have a unique lifting to E[xi, . . . , x„]p. 

Since E[a;i, . . . , Xn] is regular the Jacobian image 7r(^7) has rank I equal to the 
height of the prime. Let (ji, . . . ,ji) be the column indices of a non- vanishing I x l- 
minor and choose M := (e^j, . . . , e^j). Here Cj denotes the column vector with 1 in 
position j and otherwise. Since multiplication by M extracts exactly the columns 
of this minor also 'n{J)M has rank I. 

Applying the above formula to the components of the vector [fi, . . . , fr)^ of 
generators of p we find that its 7r-image is of the form 

(7r(/i), . . . , TT{fr)f + ^{J)M (ti, . . . , tif + (higher order terms). 

^ V ' V ' 

= (0,...,0) =:L(/i,...J,) 
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The rank condition assures that the components of L{fi, . . . , fr) generate m/tn^ as 
Fo-vector space. Hence m = (0(/i), . . . ,4>{fr)) by Nakayama's lemma. 

Now let F C E[xi, . . . , a;„]p be a coefHcient field in the sense of Theorem A.I., 
i.e., restricting the canonic projection to F gives an isomorphism 

F ^ E[a;i, . . . , a;„]p/pE[a:;i, . . . ,a;„]p = Fq. 

But then also (/)(F) is a coefficient field of Fopi, . . . Now we proceed as in 
the proof of Theorem A.l. The paragraph above shows that 4> is surjective and 
comparing dimensions one proves injectivity. □ 

To show correctness of our algorithms we need that completion commutes with 
two common operations, namely, building the integral closure and computing the 
blow up algebra. 

Lemma A. 3 (Integral Closure and Completion). Let A be a finitely generated 
E-algebra, f e A[z] a monic polynomial, p (Z A a prime ideal and consider the 
following diagram: 

IC(^[z]/(/)) -IC(Ap[z]/(/)) 

A ^i; 

Then q i-^ q' := q IC(Ap gives a bijective correspondence between prime 

ideals q C lC{A[z]/ (f)) above p and prime ideals q' C IC(Ap [z]/(/)) above pAp, and 

the induced homomorphisms lC{A[z]/ {f))q —>■ lC{Ap[z]/ {f))qi are isomorphisms. 

Proof. The fact that building the integral closure and completion commutes is 
known as Zariski's Main Theorem, see, e.g., [40, Thm. VIII. 33]. □ 

Lemma A. 4 (Blowing up and Completion). Let A be a finitely generated K-algebra, 
p C A a prime ideal and consider the diagram 

B ^ B' -.^TpdjAB 



A ^Ap 

where B is a coordinate ring of an affine chart of the blow up of Spec A at p. 

Then B' is a coordinate ring of an affine chart of the blow up of Spec Ap at mAp . 
Further q q' := c{B' gives a bijective correspondence between prime ideals q C B 
above p and prime ideals q' C B' above mAp, and the induced homomorphisms 
Bq — > B'^i are isomorphisms. 

Proof. This is an algebraic transcription of [23, Prop. A. 14. 7] which holds analo- 
gously for ground fields which are not algebraically closed. □ 

Finally we state the theorem which provides the theoretical basis for the whole 
formal desingularization procedure. 

Theorem A. 5 (Theorem of Jung-Abhyankar). Let f G E[a;i, . . . , a;„][z] be a 
monic, squarefree polynomial s.t. disCz(/) = ■ ■ ■ x'^u{xi, . . . ,x„) where m < 
n and u(0, ...,0) ^ 0. Then there is a natural number d > 1 and there are 
deg^{f) distinct power series ai € Efx^'^ , . . . , xU.'^ , x^+i, . ■ . , Xnl solving f, i.e., 
/(xi, . . .,Xn,ai) = /or 1 < i < dcg^{f). 
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Proof. An irreducible factor of / must again be monic and its discriminant must 
be a factor of discz(/), hence, can also be written as above. Now the statement 
can be found in [23, Prop. 3.2.5]. □ 
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